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1 Introduction 



Erwin Schrodinger called entanglement the characteristic trait of quantum mechanics, the 
single concept that stands out not to have any counterpart in classical physics. Richard 
Feynman said that there is only one mystery in quantum mechanics, and that is entan- 
glement and the consequences thereof. Albert Einstein was - along with many of his 
contemporaries - uncomfortable with what he called "spooky action at a distance", and 
co-wrote the seminal paper [EPR35J where the authors dismissed the idea that quantum 
mechanics could possibly be a complete description of Nature on the grounds that it con- 
tained entanglement. Today entanglement is accepted as a trait of Nature, and though 
many aspects are poorly understood, the idea of a quantum computer has shifted much of 
the research away from fundamental investigation to pragmatical work on how to utilize 
entanglement in quantum computing and quantum communication [ NCOO] . 

Entanglement is essentially a question of correlations. Correlations are a well known daily 
life phenomenon. We all know that wearing a seat belt is correlated with our chances of 
avoiding serious injury in a car accident. Other correlations are less obvious, and are often 
the focus of research, e.g. whether there is a correlation between the CO2 content in the 
atmosphere and global warming, whether there is a correlation between certain diets and 
the risk of acquiring cancer, or how correlated the number of years of education is with 
income. In classical theories a correlation can never be better than perfect, the best we 
can do is a perfect correlation, e.g. if a certain number of years in school linearly increases 
the yearly income. As all scientists know, this is not the case, but probably there is a less 
than perfect, though positive correlation. However, quantum mechanics opens the case 
that there might be better than perfect correlations. This is the case when considering e.g. 
the spin of two electrons separated by some spatial distance. We can measure the spin of 
one, and it may not only be correlated with what we measure on the other electron, but it 
will be more than 100% correlated; the electrons were entangled before the measurement. 
In the simplistic case described here, this will never be evident, all we can see in our 
measurements is a perfect correlation. However, John Bell wrote down an inequality in 
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Figure 1.1: A physical implementation of "Bob", one of the characters of quantum information 
theory. This is the setup used in an optical quantum cryptography experiment at the Department 
of Electronics and Telecommunications, NTNU. Image courtesy of Vadim Makarov, used with 
permission. 

1964 |Bel64] that would be satisfied for all classical correlations, but would possibly be 
broken if there were quantum correlations present. Now it is accepted as an experimental 
fact that Bell's inequality may be broken, as first confidently demonstrated by Alain Aspect 
et al. in 1982 [ADR82] . Thus Nature seems to allow quantum correlations on top of the 
classical correlations. This excess correlation is what we know as entanglement, and is 
evident whenever Bell's inequality is broken. 

Though the connection is not obvious, entanglement is today mainly studied in connection 
with quantum information theory, as it is believed that a quantum computer will - if ever 
built - utilize entanglement to perform certain computational tasks faster than classically 
possible. Feynman pointed out in 1982 that simulating quantum systems on a classical 
computer was nearly impossible due to the much larger Hilbert space of the quantum 
variant |Fey82| . Thus, he suggested the use of a quantum computer to simulate quantum 
systems. However, it was not until Peter Shor's discovery in 1994 |Sho94] that a quantum 
computer could factorize large numbers exponentially faster than classically possible that 
the idea really gained momentum. The quantum computer is still very much a theoretical 
device, despite considerable experimental effort on small scale systems. It has so far not 
emerged any single clear candidate for the implementation of the quantum bits (qubits) 
and the operations needed to process the quantum information. The NMR experiment 
that factorized the number 15 in 2001 !\'SB '01 was the first implementation where a 
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quantum algorithm was proved to work, though on a relatively trivial problem. The NMR 
implementation, along with most implementations based on atomic, molecular or optical 
physics (AMO), is poorly scalable. The approach based on ion traps has nevertheless 
become fashionable, and was one of the first promising proposals for a quantum computer 
made in 1995 by Cirac and Zoller [CZ95J. Other implementations are based on condensed 
matter systems, such as superconducting qubits or quantum dots. The main problem in 
many of these approaches is decoherence, that is the qubits become entangled with the 
environment and thus information is lost into the environment during the computation. 
Recently, two superconducting qubits were entangled [SAB + 06] , making it an ever more 
promising candidate for quantum computing. Also recently, a hybrid approach has been 
proposed, exploiting the advantages of different schemes [ADD + 06j . where trapped polar 
molecules are controlled by a microwave field set up by running superconductors. Despite 
the extensive effort, one still doesn't know whether the quantum computer even can be 
built, or if one might encounter a fundamental reason why this is an impossible feat. Even 
so, the process will learn us fundamental facts about physics and the process cannot be 
said to have been in vain. 

Phase transitions are also familiar phenomena of daily life, like the freezing of water at 
zero temperature. Such classical phase transitions occur when the temperature is changed 
across some critical value T c , and the thermal fluctuations are exactly what drives the 
system from one thermodynamical state to another. Quantum phase transitions on the 
contrary, can happen at zero temperature, and therefore without thermal fluctuations. 
These continuous phase transitions occur when the ground state changes nature through 
the variation of a parameter, such as an external magnetic field. The only fluctuations 
present are the quantum fluctuations that must exist due to the Heisenberg uncertainty 
relation, and hence the name quantum phase transitions. The wave function of the system 
at the critical point is very complicated, as is evident from the fact that it carries more 
entanglement that at off-critical points. The complexity of the wave function at the critical 
point also means that more information is encoded in the wave function, thus making it 
more usable for quantum information. 

This thesis is organized as follows; Chapter [2] is a brief introduction to the most important 
principles in entanglement theory, with an emphasis on those subjects which are considered 
in later chapters. Chapter [3] addresses some general issues when it comes to modelling 
quantum systems on a classical computer, and the line is drawn between renormalization 
schemes and the problem of finding the ground state of a quantum system as addressed 
in |Skr06a| . In Chapter H] we investigate bosonic systems, and compute the entanglement 
therein. In particular, the results from |Skr05aj are connected to other results from the 
theory of entanglement in continuous variable systems. Finally, Chapter [5] contains the 
bulk part of the thesis, summarizes quantum critical systems, and applies ideas from |SQ05l 
ISkr05bl ISkrOfml ISkr06b| . 



In all chapters, unless otherwise specified, units such that H — c — 1 are used, where c is 
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the speed of light in vacuum. All notation is as far as possible applied as is conventional 
in the literature. In particular, the symbol © is used for the direct sum of matrices, e.g. 
A © B = diag(A,_B), while © designates direct product, which for matrices is defined in 
e.g. [NCQOj . 



2 Entanglement 



In its simplest form, entanglement arises from the possibility opened by quantum theory 
to write down a wave function that is non-local, encompassing two constituents separated 
by arbitrarily large distances. Consider a wave function of two particles shared by the 
protagonists of quantum information theory, Alice and Bob. Now, any measurement on 
either's part implies the collapse of the wave function at the other's side. This collapse 
happens instantly (in some frame of reference), and this is the source to its counter- 
intuitivity, not to mention the apparent contradiction with special relativity. However, 
Alice and Bob can never exploit this effect to communicate any information faster than 
the speed of light. Due to the statistical nature of the outcomes, the parties need to 
communicate classically to compare their results. Moreover, the simple scheme here can 
easily be explained through classical means with some "hidden variable" interpretation. 
However, there are more complicated schemes that efficiently demonstrate the breaking 
of Bell- type inequalities, thus ruling out any hidden variable explanation [Mer85j. In this 
thesis we will restrict our investigation to bipartite entanglement, and ignore the case where 
there are more parties than Alice and Bob. Thus, consider two systems in separate Hilbert 
spaces TCa and TCb respectively. Given a state vector \i/))a £ 7~La in the first system and 
a state vector \<p)b £ "Hb in the second, we can define a state in the entire Hilbert space, 
H = Ha <8> Hb as 

Inseparable = \lf>) A ® \<P)b- 

As indicated, this is a separable state, and any state that can be written in this way is 
separable, and thus carries no entanglement. More generally, the state of the entire system 
can be written ClS db Schmidt decomposition 

x 

\^)=^V^i\A)A®Wi)B (2.1) 

8=1 

where the upper limit of the sum, x? is called the Schmidt number of the state, and 
is maximally equal to the smallest of the dimensionalities of the two Hilbert spaces in 
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question. The states are orthogonal for all i, and the same is true for \(pi). The 
expansion is unique, and the Schmidt number can be considered a very brute measure of 
entanglement, in particular the measure E(\^f)) = log 2 x fulfills some of the requirements 
for an entanglement measure as described later |Vid03| . 

Identifying an entangled state can, as opposed to what impression one gets from the simple 
analysis above, be a very hard task. The problem can be stated as finding the basis in 
which one can write a given state as a separable one, or proving that no such basis exists. 
It is important to realise that this task depends heavily on the choice of partitioning of 
the system. The task of writing down this basis indeed seems so daunting that one often 
resorts to simply stating whether a state is separable or not, and if it is non-separable, how 
much entanglement do we have? We will return to this problem in Section I2TT1 

Given a Schmidt decomposition of a pure state \^f), one can straightforwardly write down 
the density matrix of (say) TCa as 

i 

where we have omitted the index A since the identity of the Hilbert space is obvious. This 
density matrix contains all information accessible to Alice. Hence, finding the Schmidt 
decomposition of the pure state wave function is equivalent to diagonalizing the reduced 
density matrices of the two subspaces. That is one indication of the importance of the 
Schmidt decomposition. 

More generally though, the state in the full Hilbert space 7i may not be a pure one, but 
an indeterminate mixed state with density matrix p. Still, the reduced density matrix 
of H. a is well defined, but the entanglement is not, in the sense that one cannot find a 
genuine measure of entanglement for the state. Purification learns us that any mixed state 
in a Hilbert space 7i can always be interpreted as a pure state in some (imaginary) larger 
Hilbert space H' D H. Consequently, the physical state p G H is entangled to the system 
of the larger Hilbert space H'. 

If Alice has access to the pure states |p!^ j e H-a and Bob has the pure states j/4 J e 
Hb, the state 

P = £ P*J ® pf ] Yl Vi i = 1 

ij ij 

carries only classical correlations, and must thus be classified as unentangled, or separable 
|Wer89j . A density matrix that is non-separable consequently is entangled, but finding a 
decomposition for a given density matrix is computationally hard. Moreover, for the case 
where p is a mixed state, it is not yet known how to measure the entanglement of the state 
such that classical correlations are eliminated. 

Assume Alice and Bob being at different locations, separated by a large distance, and shar- 
ing a state describe by the density matrix p. Now, each party may operate on their state 
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by unitary transformations or measurements, and they are even allowed to communicate 
classically, e.g. through telephone, but any joint operation on both parties are prohib- 
ited. This protocol is known as LOCC, Local Operations and Classical Communication. 
Through LOCC one can generate all classical correlations, but no entanglement can ever 
be created this way. This means that LOCC describes the fundamental difference between 
what can be created classically and quantum mechanically and is a useful starting point to 
describe quantum correlations. What is sometimes called the fundamental law of quantum 
information processing can be phrased as 



Through LOCC alone, Alice and Bob cannot increase the total 
entanglement of the state they share. 



Mathematically, local operations can be either unitary operators Ua G *Ha and [/gS Tis 
or general measurements |a| an d {^}" ^ e g eneranze d measurements are complete, 
in the sense that 

i i 

where Aj are measurements performed exclusively on Alice's system, and Bi are performed 
on Bob's. After Alice has applied a measurement on her state p^ A \ihe state collapses into 
p' = Aip^ A\/p(i) with probability p(i) = Tr (Aip^A\), and correspondingly for Bob's 
measurements. LOCCs are important in the sense that it should be assumed that entangle- 
ment cannot increase under LOCC, and should be invariant if there are no measurements 
involved. 

Moreover, a fundamental fact involving LOCC is that, given two states j^i) and |^ 2 ), of 
the joint system of Alice and Bob, one can transform the first into the latter by LOCC only 
if the ordered eigenvalues Ai of p^ = Tr is majorized by the eigenvalues A 2 of 

^(A) _ rp r B |^ 2 ^^ 2 |. A decreasingly ordered vector v majorizes another ordered sequence 
w, written v >- w if 

k k 

Vi > Wi for all k — 1, • • • , d 

i=l i=l 

where d is the smallest of the dimensions of the two vectors. 



2.1 Entanglement measures 

As we have pointed out, entanglement is very much considered a fundamental resource in 
nature, much like energy or information are considered resources in other contexts. And 
anything called a resource must be measured, so we would have some idea about how 



8 Entanglement and Its Applications in Systems with Many Degrees of Freedom 



"much" we have of that resource. Thus one needs to define entanglement measures, a very 
non-trivial task Myh04|. For some state described by the density matrix p, we want the 
entanglement E(p) G R of the state. This measure must fulfill some basic requirements to 
be physically reasonable |VPRK97l iHHHOOj ; 

El E(p) > with equality if p is separable. 

E2 E(p) = E(U A <g> U B pU ] A <g> U B ) where U A and U B are any local unitary operators on 
Hilbert spaces Ha and TLb respectively. 

E3 The entanglement cannot increase under LOCC, E(<dp) < E{p) where G is any 
LOCC operator. 

E4 Convexity: Y,iPi E (Pi) > E (I2iPiPi)- 
E5 Partial additivity: E(p® n ) = nE(p). 

E6 Continuity: If lim rwoo (?/>® n |p n |'i/> (g,ri ) = 1 with p n a density matrix of n pairs, then 

lim-[£(|^r)-£(Pn)]=0 

n— >oo 77, 

One also often requires the normalization that the Bell state \(3) = ( 1 00) + |11)) / V2 should 
have unit entanglement. Any entanglement measure that satisfies conditions El-4 is known 
as an entanglement monotone. 

2.1.1 Forming and distilling entanglement 

Given some quantum state p, the question of how much entanglement this contains is 
nontrivial. However, conceptually the question can be simplified by defining entanglement 
of formation and -distillation. Assume that we have a large number of Bell states, and are 
allowed to use only LOCC on these states, how many copies of a given bipartite state p 
can we make? The number of copies of p divided by the number of Bell states is known as 
the entanglement of formation Ep when one takes the limit of an infinite number of Bell 
states. On the contrary, the entanglement of distillation is defined as the number of Bell 
states we can create from a large number of copies of p. The distillable entanglement Ep, 
is always less than or equal to the entanglement of formation since one inevitably loses 
information in the process of formation and distillation. Moreover, it turns out that any 
reasonable entanglement measure is bounded by these, [HHH00J 

E D {p) <E(p)< E F (p). (2.2) 

For the special case of p = |^)(^| being a pure state, it turns out that one can uniquely de- 
fine a bipartite entanglement measure since the entanglement of formation and -distillation 
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are equal and equal to the von Neumann entropy of the reduced density matrix of one of 
the parties. That is, 



EfMM) = S(p') = -Trp'log 2 p' with p' = TirslVWI 



(2.3) 



where p' is the reduced density matrix accessible to Alice. The entanglement must be 
symmetric, so it does not matter whether we trace out Alice's or Bob's system. The von 
Neumann entropy is nothing but the Shannon entropy of the eigenvalues of the reduced 
density matrix, as the Shannon entropy of some random variable X with probabilities 
a:i,...,a;jv is H(X) = —J2 n x n log 2 x n . Note that we use logarithm base two, meaning 
that the von Neumann entropy is measured in ebits (entanglement bits), and the von 
Neumann entropy of a shared Bell state is unity. 



Consider the von Neumann entropy in the case of a single qubit, that is the density matrix 

p = a \0)(0\ + 010X11 + F\1)(0\ + (1 - a)|l)(l| 
This has two eigenvalues A and 1 — A, and thus the von Neumann entropy becomes 



a (3 
[3* I- a 



S{p) = -A log 2 A - (1 - A) log 2 (l - A) = H(X), 



(2.4) 



which is known as the binary entropy function, shown in Figure 12.11 H(X) is zero when 
A = 0, 1 and reaches maximum at A = 1/2 as H(l/2) = 1. Obviously, the former is a pure 
state, while the latter is a maximally entangled state for a single qubit. 

For a general mixed state, the situation is more complicated. Indeed, it is very hard to 
eliminate the classical correlations in a mixed state such that the entanglement measure 
contains only the quantum correlations. The entanglement of formation for a mixed state 




Figure 2.1: The binary entropy function H(\) = — Alog 2 A — (1 — A)log 2 (l — A), which is the 
entanglement of a single qubit whose density matrix has eigenvalues A and 1 — A. 
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p is defined as 

E F (p)=milj2^ S (\^)(M |p = ^A fe |^)(^| I (2.5) 

I k k ) 

where the infimum is to be taken over all decompositions of p into pure states. Here, 
S( |) indicates the von Neumann entropy of (say) Alice's reduced density matrix 
p' = Tr B\ip)(ip\. The explicit computation of the entanglement of formation is known only 
for a limited number of cases, such as for the two-spin entanglement known as concurrence 
reviewed next, and the two-mode Gaussian entanglement reviewed in section | 



2.1.2 Concurrence 



Concurrence is a measure of entanglement in a system consisting of two spins, or any other 
two- level system that is closely related to the entanglement of formation |HW97} IWoo98] . 
We consider a mixed state p, such as the reduced density matrix of two spins in a larger 
spin chain. The spin-flip operator on a qubit is defined as = <J v \ip*)., which flips the 
spin in the standard basis. Note that complex conjugation, like transposition, is a basis- 
dependent transformation and it is thus not a physical operation. Now, to flip both spins 
in the standard basis {| ft), 1 1 | [[)} take the transformed density matrix 

~p= [p v ®a v )p*(a y ®a y ), 

and define the Hermitian matrix R = \J p 1 l 2 pp y l 2 whose decreasingly ordered eigenvalues 
we denote Ai > A2 > A3 > A4, and assume that at least two of them are non-zero. The 
trace of R is simply the fidelity between the spin-flipped density matrix and the original^]. 
Now the concurrence is defined as 

C(p) = max(0, Ai - A 2 - A 3 - A 4 ). (2.6) 

The concurrence is a valid entanglement measure in its own right, but the entanglement 
of formation equals 

1 



E F (p) = Hl-[l + y/l-C(p)*\), (2.7) 

where H(x) is the binary entropy function (12.41) . Since Ep{p) is a monotonically increasing 
function of C(p) and they coincide at zero and unity, they can be considered equivalent 
entanglement measures, and one usually call the concurrence simply the entanglement of 
formation due to this equivalence. 

A conceptual sketch of the different ways to measure entanglement by concurrence and 
entanglement entropy is shown in Figure 12.21 The key point is that the entropy measures 



1 The fidelity is a distance measure between two quantum states, being unity if they are equal, and less 
than one, though positive, else. The fidelity is revisited in Section I5H1 
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the entanglement of a set of (not necessarily contiguous) spins with the rest of the chain, 
where the whole chain is in a pure state. The concurrence, however measures the entangle- 
ment shared between only two spins, where these two spins together might be in a mixed 
state. 



p> » S(p') 




p> » C(p') 

Figure 2.2: Conceptual sketch of the two entanglement measures entanglement entropy and 
concurrence in a spin chain. The former amounts to finding the reduced density matrix of a 
part of the spin chain and hence the entropy S(p) = — Tr p\og 2 p, while the second considers two 
individual spins in the chain and their entanglement is measured by the concurrence. 



2.1.3 Other entanglement measures 

There is a number of other measures of entanglement for mixed states that are entangle- 
ment monotones, most of which reduce to the entanglement entropy when one considers 
pure states. These include negativity, which is defined on a 2-qubit state as M{p) = 
2max(0, — A neg ) where A neg is the sum of the negative eigenvalues of p Ts where one takes 
the partial transpose of Bob's part of the density matrix |VW02j . 

Another measure that has gained attention recently is squashed entanglement or conditional 
mutual information (CMI) [CW04J, which is defined as 

E C mi(pab) = \ inf S(A : B\C) 

Z PABC 

where the infimum is to be taken over all tripartite states that are such that pab = 
Tr cPabc- The entropy S(A : B\C) is known classically as conditional mutual information 
and is the mutual information of A and B given C. Mutual information is given by 

S(A:B) = S(p A )+S(p B )-S(p AB ), 

while conditional mutual information is 



S(A : B\C) = S(p c ) - S(p A c) - S(p B c) + S{ PABC ). 
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That is, we introduce an artificial system C, and compute the maximal information content 
in A and B given this system. Squashed entanglement is known to fulfill many of the 
requirements for entanglement measures and is a promising candidate for the final measure 
of entanglement for mixed states. However, the definition includes an extremalization that 
may be computationally very hard, though little has been done on computing squashed 
entanglement [Chr06] . 

2.1.4 Renyi entropy 

The entanglement entropy (12.31) is not the only entropy one can define given the density 
matrix p, whose eigenvalues are \. In particular, the von Neumann entropy is a special 
case of the Renyi entropy, 



since S(p) = liniQ,^! H a (p). The Renyi entropy in a classical context where the eigenvalue 
spectrum is replaced by a probability distribution, has applications in several areas of 
physics and information theory, e.g. in signal processing [Jen05] . along with the Shannon 
entropy, which is the classical equivalent of the von Neumann entropy. All Renyi entropies 
are zero if and only if p is a pure state and positive otherwise. Finally, in the limit a — > oo 
the Renyi entropy is simply the logarithm of the largest eigenvalue, 



This is known as the single-copy entanglement [EC05J, since it measures how much entan- 
glement can be extracted from a single specimen of p, rather than the asymptotic notion 
used to define entanglement of formation that leads to the von Neumann entropy. 



Entanglement is an important subject in its own right, as it challenges our understanding 
of Nature and how it works, and might be at the heart of physics. Nevertheless, it is in 
quantum information theory that applications of entanglement are known, and hence where 
the subject has gained the most attention. Shor's algorithm is an example of a quantum 
algorithm that utilizes entanglement as a computational resource. However, entanglement 
is not necessary for quantum computation, e.g. Grover's search algorithm does not attain 
its power from entanglement [Llo99] . This might be the reason that Grover's algorithm 
is a mere quadratic speedup, compared to Shor's exponential speedup relative to known 
classical algorithms. 





Hoo(p) = -log 2 maxAi. 



(2.9) 



2.2 Entanglement in quantum information 
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Also, entanglement can be used as a resource in quantum communication, including super- 
dense coding, where two classical bits can be transmitted by one qubit given that Alice and 
Bob share an entangled state before the communication |BW92] . Quantum teleportation 
[BBC + 93j exploits a shared entangled state to transmit an arbitrary quantum state from 
Alice to Bob using only LOCC, a technique that has been experimentally implemented 
in atomic systems to much popular interest |KE04j . The quantum cryptography scheme 
due to Ekert also relies on distribution of entangled particles jEke91] , though the popular 
BB84 protocol does not use entanglement [B B84] . 



3 Numerical Analysis 



The study of quantum mechanical systems that involve more than a few particles inevitably 
requires the use of classical computers to solve the system. The size of the Hilbert space of 
a system of N qubits increases exponentially as 2^, which is a contributing factor that a 
quantum computer may attain its vast computational power, but also naturally hampers 
simulations of quantum systems on a classical computer. And as long as we don't have 
the quantum computer, the simulations must be done on the classical variant. Given a 
quantum system of TV spin-1/2 particles, a possible encoding of the qubit, a separable state 
is 

N 

l^sep) = ®(a n |0> n + /? n |l) n ) (3.1) 

71=1 

with a n and f3 n complex coefficients. This gives 2N degrees of freedom — overall phase 
factors and normalization corrected for. However, a general state is 

2^-1 
r?=0 

where \rj) are the orthogonal basis states, and ^ \a v \ 2 = 1. The as now constitute 2^ — 2 
degrees of freedom. Hence, when entanglement is not included, the complexity in the state 
is dramatically reduced as opposed to the general entangled state, although the schematic 
overview here is somewhat simplified. 

It is also immediately clear that if one computes the spin-spin correlation function 

for the separable state (13. ip . it is identically zero. Hence, if there are any correlations in 
the state, they must arise from entanglement unless there are degeneracies or finite tem- 
perature, which means we may not have a pure state. It is a well-known fact from critical 
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phenomena that the correlation function decays exponentially with inter-spin distance, 
Tij ~ for large distances, \i — j\ ^> 1. £ is then known as the correlation length. 

However, at critical points, the correlation function decays polynomially, thus demonstrat- 
ing the fact that entanglement is much more prevalent at critical points than elsewhere in 
the parameter space. 

Sometimes it is possible to map the system of interacting qubits onto a simpler system 
where the degrees of freedom grows linearly as opposed to exponentially, as we will explore 
in later chapters. However, this is not always possible, and much attention has been drawn 
to the subject of finding numerical recipes that can be used to efficiently compute e.g. the 
ground state of a Hamiltonian. 



3.1 Using symmetries 

Most models exhibit symmetries, that is operators that commute with the Hamiltonian, 
and thus share a common eigenvalue set with this. Imposing periodic boundary conditions 
on a model will lead to translational invariance, and thus conservation of momentum 
according to Noether's theorem. Thus, the full Hilbert space of the Hamiltonian can be 
reduced into subspaces with definite momenta, and thereby somewhat reducing the size of 
the matrices that one needs to diagonalize. This is not as efficient as other methods - 
by a long shot — but it may reduce the complexity enough to allow numerical analysis of 
slightly larger systems than otherwise possible. Nevertheless, the size of the matrices still 
grows exponentially and therefore no fundamentally new window of exploration is opened. 

We focus on fermionic models, that is where a Fock state wave function can be written as 
a string of binary digits, and we consider translational, reflection, and parity invariance. 
Define the operators T, TZ, and V for these operations respectively. Each of these symmetry 
operators S e {T, 1Z, V} commute with the Hamiltonian, [H, S] = 0, and each has an order 
n for which S n = 1. Thus, S must have n eigenvalues = e 2nlk / n _ The symmetry operators 
and their eigenvalues are summarized in table 13.11 Each of these operators will reduce the 
Hilbert space into their respective subspaces, thus reducing the dimensionality accordingly. 
However, they do not commute with each other, e.g. [T, TZ] ^ 0, so it may not be possible 
to reduce the complexity with all operators simultaneously. However, [P, T] = 0, so we 

Table 3.1: The symmetry operators of a generic spin chain and their eigenvalues and -states. 
\x) is an eigenstate of the Hamiltonian. 



Name 


Symbol 


Order 


Eigenvalues 


Eigenstates 


Translation 


r 


N 


rjk = exp(2iiik/N) 


En^o e 27rikn/N T n \x) 


Reflection 


K 


2 


r = ±1 


(l±W)|x> 


Parity change 


V 


2 


p = ±1 


(i±P)lx> 
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exploit these two symmetries simultaneously to block-diagonalize the Hamiltonian. The 
size of the matrix for a given partition of the Hilbert space with quantum numbers (k,p) 
is denoted N/ k;P y 

If we look at the translation operator, this will provide the greatest reduction of the Hilbert 
space since it has the largest number of different eigenvalues, while the parity operator will 
split the Hilbert space in two. Given a state \x), not all values of k are allowed though, 
e.g. look at the state \x) = 1 0101) ; 

(1 + e 2 - fc / 4 r + e 4 - fc / 4 r 2 + e 6 - ifc / 4 r 3 ) ix) 

= (1 + e" ifc ) |0101) + (e" ifc/2 + e 3 " ife/2 ) |1010) 
= (1 + e" ifc ) (|0101) + e" ifc/2 |1010)) . 

The latter expression is zero unless k = or k = 2, and hence k = 1 and k = 3 do not 
correspond to eigenvalues of the operator. More generally, if T s \x) = \x) f° r an s < iV 
and some state \x), then N/s must be an integer, and k must be a multiple of this integer 
to be a valid quantum number. The two extreme cases are if s — N, e.g. |0001), when 
all k = 0, • • • , N — 1 are good quantum numbers and if s = 1, e.g. |0000), when only 
k = is a good quantum number. Hence, for a given state \xi), there exists s x valid 
fc-quantum numbers, where s x is the number of unique states that \x) can be transformed 
into through the translation operator. We have observed that the ground state always is in 
the k = segment of the Hilbert space, at least assuming there is no spontaneous violation 
of translation invariance. 

To exemplify, the N = 4 Hamiltonian is a 16 x 16 matrix. This can be reduced according 
to Table 13.21 into groups with a defined parity and wave number. As we see, the matrix 
is reduced into one 4x4 and smaller matrices. In any one of these eight groups, each 
generating state \x) defines a translationally invariant state, 

|^(M) = e 2wikn/N T n \ X ). 

Now each element of the new Nn.^ x Nn.^ Hamiltonian matrix h is 

s x — 1 s x '-! 

/»« = (^\H\^) =J2Y, (x\T n HT m \x')^ lk{n - m)IN 

n=0 m=0 

for two generating states |%) and |%'). Hence the computational complexity is reduced 
somewhat. Figure 13.11 shows how large the submatrices are for different system sizes. 
The figure shows that with a computational threshold that restrict matrix diagonalization 
upwards to A^ = 10 4 , the size of the system that can be considered is increased from A^ = 13 
to N = 18, which can be a crucial improvement in certain context, such as the criticality 
detection discussed in Section 15.4.51 Also, one may go beyond this threshold when looking 
at individual points in the parameter space, but if scanning through a large portion of the 
space one would usually apply smaller system sizes to avoid excessive time consumption. 
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Table 3.2: The groups when using translational and parity symmetries in a fermionic spin chain 
of size N = 4. 



Quantum numbers 


X) s x N (k, P ) 


k = 0,p=l 


0000) 1 
0011) 4 
0101) 2 

mi) l 


k = 1, p = 1 


0011) 4 1 


k = 2, p= 1 


0011) 4 
0101) 2 


k — 3, p — 1 


0011) 4 1 


fc = 0, p= -1 


0001) 4 
0111) 4 


fc = 1, p = -1 


0001) 4 
0111) 4 


fc = 2, p= -1 


0001) 4 
0111) 4 


fc = 3, p = — 1 


0001) 4 
0111) 4 




Figure 3.1: The reduction in complexity when utilizing symmetry as measured by the size of 
the matrices needed to diagonalize. The dashed line is 2^, which is the size of the matrix if no 
symmetries are used, while the points are the size of the matrices in each (k,p) segment of the 
Hilbert space. For illustration the horizontal line N( k ^ = 10 4 is shown as a threshold for how 
large matrices we can diagonalize, and the improvement when using symmetries is from N = 13 
to N = 18. 
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3.2 Renormalization techniques 

Renormalization in classical condensed matter systems was conceived by Leo P. Kadanoff in 
1966 [Kad66j, while Wilson |Wil75j refined the method and applied it to critical phenomena 
such as the Kondo effect. The basic realization is that at critical points, the system becomes 
scale invariant, that is the thermodynamic properties can be derived from what happens 
on large length scales. Hence, one can construct a mapping of a group of lattice sites onto a 
"supersite" whose properties are derived from the smaller group. This can now be iterated 
to successively larger blocks until the thermodynamic properties can be extracted. This 
is known as the renormalization group (RG), and it creates a renormalization flow in the 
parameter space. The technique has proved highly successful in classical lattice models, 
and has become a part of the standard theory of critical phenomena. 

However, the RG method fails spectacularly for some very simple models, such as the one 
dimensional particle in a box problem. The very problem that is often used as the first 
example in quantum mechanics textbooks could not be solved by this approach. White 
and Noack |WN92] identified the crucial point, that the ground state wave function which 
is a sine wave, cannot be composed of the ground states of smaller partitions which are 
also sine waves but with different end points. White thus proposed a new renormalization 
group, the density matrix renormalization group (DMRG), which picks out the relevant 
wave functions to keep in further iterations very effectively [Whi92j . In brief, the technique 
amounts to dividing the full system into a "system" , an "environment" , and two single sites 
between them that work as a connection. Now, the ground state of this so-called superblock 
is found as a first approximation. The procedure is then iterated by shrinking the system 
until it consists of a single site and then growing the system until it reaches a threshold 
size, and so on until the desired accuracy is reached |Sch05] . This has proven successful 
for a number of models, such as the above particle in a box and more real problems such 
as the Heisenberg model, and has established itself as the predominant method for solving 
low-level excitations in one dimensional systems. 

However, the DMRG fails to reproduce the algebraic decay of the correlation functions 
at the critical point, which is a crucial requirement for any such technique. It has been 
conjectured that this is because, while the DMRG technique takes some of the entangle- 
ment into consideration, it does not preserve maximal entanglement under renormalization 
[ON02bj . At a critical point it seems that the entanglement constitutes a major part of the 
wave function, such that an RG technique to estimate e.g. the ground state wave function 
here must take entanglement into account. 

There have been a number of suggestions surfacing in the latest years, such as the en- 
tanglement renormalization [V7d 05j . which tries to devise a new renormalization scheme 
based on projecting out the most entangled states. Another approach based on so-called 
projected entangled pair states (PEPS) |VG04l IVWPGC06j enables efficient computation 
in two or more dimensions. 
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Some issues concerning entanglement properties of the ground state in the Ising chain are 
addressed in [Skr06aJ . In particular, we see that even though the entanglement is a highly 
important attribute of the wave function, and especially so at the critical point in the 
Ising chain, there are very few terms in the Schmidt decomposition that needs to be taken 
into account to construct the wave function. However, it is not yet clear how one should 
proceed to identify these terms a priori, that is when the wave function is not previously 
known as is the case with the Ising chain. Nevertheless, it is enough to know the first four 
terms in the Schmidt decomposition of a spin chain with iV = 100 spins to determine the 
entropy with an accuracy of 10 -4 . 



4 Continuous Variables 



Usually quantum information theory is occupied with two-level quantum systems, both due 
to its simplicity and its obvious connection to classical information theory. In particular, 
the analogy that classical bits are substituted with quantum bits in a quantum computer 
is appealing, but quantum information in systems with more than two levels, or even a 
continuous spectrum, is also considered. Such systems are experimentally important since 
they are easy to generate and manipulate, and quantum information based on quantum 
optics is arguably where most effort is done implementing feasible quantum systems. Par- 
ticularly, optical realizations are important in quantum communication. In this chapter we 
will study continuous variable (CV) systems, that is where each entity can have an infinite 
number of eigenvalues rather than the two-level system of a qubit. The infinite number 
of eigenvalues corresponds to the fact that any given mode of the system may contain an 
arbitrary number of particles, denoted \n) with n — 0, 1, 2, • • • . The canonical example of a 
CV system is the textbook harmonic oscillator, which has an infinite set of discrete energy 
eigenvalues E n = u(n + 1/2). The density matrix of such a system is, in contrast to the 
fermionic case, infinite dimensional. This means a very different approach to computing 
physical properties, and also that the entanglement of a state is not bounded upwards as 
opposed to the fermionic case. 

Generally, a CV system of N particles, or modes, consists of iV conjugate pairs x k and p k 
that satisfy the canonical commutation relations 



[x k ,pi] = iSki. 



(4.1) 



It is customary to define the annihilation and creation operators 





(4.2) 



which satisfy the commutation relations [a n , a, 
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The conventional way to define the density operator is in terms of x and p, p(x,p). We 
can also explicitly write down the density matrix in the position basis, in which it is 
a continuous matrix p(x,x') = t\)* n (x)p(x', p')i/j n (x') in a choice of basis with basis 
vectors {ip n (x)}. This normalizes as J p(x,x) d 2N x = 1 and the expectation value of an 
operator 0(x,p) is 

(0(x,p)) = TrO(x,p)p(x,x'). (4.3) 

However, it is not always convenient to write the states in the explicit position represen- 
tation, and we will review the equivalent Weyl representation of the states below. 

Separability criteria for CV systems have been proven [GKLCOlJ, and we will here recite 
some of the important aspects of the continuous variable systems in terms of Gaussian 
states and refer to the proofs in the literature. 



4.1 Weyl algebra 

Continuous variable states as defined by the Xk and pk operators can be conveniently 
formulated in terms of the Weyl algebra. To this end, we define the 2N operators as 

r = [xx, pi, x 2 , P2, ■■■ x N , p N ] T . (4.4) 
These operators' commutation relations can be written as 

[r k ,ri\ = Um (4.5) 

where J is the 2N x 2N matrix 

N 



k=l 



1 

-1 



(4.6) 



This matrix is the defining matrix of the symplectic group Sp(2iV), which is the group 
under multiplication of all 2N x 2N matrices M that satisfy 

M T JM = J. (4.7) 
Now we may define the Weyl operators which are defined on £ £ R 2JV as 



W(0 = expf-i^£ fe rfcj 



(4. 



Finally, defining the form cr(£, £) = ^iJijCji the Weyl relations are 

MOMC) = e- iCT ^w(c)w(0. ( 4 - 9 ) 
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Given a quantum state p, the characteristic function x{0 °f the state is simply the expec- 
tation value of the Weyl operators, 

= Tr P W(0. (4.10) 

Inversely, the density matrix can be written in terms of the characteristic function and the 
Weyl operators, 

P = tAf f x(0M"0 d w £. (4.11) 



(27T) 



All expectation values of polynomials of the canonical operators can be expressed in 
terms of differentiations of the characteristic functions. Indeed, the two first momenta are 

(r k ) = i-^x(th)\ t=0 (4.12a) 
d 2 i 

^ = ~dhdf 2 e ~ 5 + t 2 e l ) , (4.12b) 



tl=*2=0 



where is the unit vector of element k. 



According to the above, the formulation of the theory with Weyl operators is equivalent 
to the formulation with canonical operators. Hence, which one is applied is a matter of 
taste and computational convenience. A general state in the form of a density matrix is 
completely determined by all nth order momenta. However, we will focus on the impor- 
tant class of Gaussian states. These are experimentally important due to the simplicity 
involved in creating and manipulating them. Also, they are simple and elegant to compute 
analytically but nevertheless involves a rich structure. Further, for the Gaussian states 
only the first and second order momenta are necessary to completely determine the state. 
This is of course a convenient situation where both experimental and theoretical simplicity 
meets and Gaussian states are therefore by far the most studied continuous variable states 
in both camps. 



4.2 Gaussian states 



A state p is called Gaussian if its characteristic function is a Gaussian, 



X(0 = exp 



(4.13) 



The 2N x 2N matrix T must be positive and symmetric to give rise to finite expectation 
values. T is known as the correlation matrix and 5 G R 2Ar as the displacement vector. 
Equivalently, and according to the convention used in |Skr05a] . a Gaussian state can be 
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written in the position basis as 



det (A' - C) 



7T 



D 



cxp 



< Mr r 'i) V, 



x exp 



- (f/,-. !,-//; + q'iA^q'j) + qiQjq'j + + d*q[ 



(4.14) 



Here A' = Re A, C = ReC, and d! = Red with A and C N x N matrices and d an iV 
dimensional vector. The equivalence of the two formulations is evident when we compute 
the first and second order momenta from Eqs. (I4.12p and from direct calculation on the 
density matrix. To simplify the equations we consider the variational matrices 



Qm = (qkqi) - (qk)(qi) (4.15a) 

Pki = (PkPi) - (Pk)(pi) (4.15b) 

Sm = ^(qkPi +piqk) - (qk)(pi)- (4.15c) 

Computing these for both formulations, we find 

Qki = T 2k ^ 2l ^= l - {(ReA-ReCy 1 } kl (4.16a) 

Pki = T 2K2l = A kl - (A kt - C^Q^Afl - C j7 ) T (4.16b) 

Su = ^r 2 fc-i, 2 * = - Qki(I™.An + ImCy (4.16c) 

(g*> = 6*-i =2Q kl Redi (4.16d) 

(Pk) = 6fc = - 2(ImA fci - \mC ki )Qij Redj + lmd k . (4.16e) 



This shows the mapping between the two formulations, and though they are equivalent, 
the mapping is complicated. However, the position representation gives a more physical 
view that counters the mathematical convenience of the Weyl representation. 

As shown in [Skr05aJ, using the position representation one can perform a rotation, scaling 
and another rotation of the coordinate system to separate the density matrix into single 
particle density matrices for virtual particles. Having found this transformation, it is 
easy to find the eigenvalues of the single particle density matrices and thus the eigenvalues 
of the full density matrix. Also, we show that the massless Klein-Gordon field possesses 
conformal symmetries (see Section IB~2l) . and that this can be identified by the entanglement 
signatures in the model. 

It is easy to realize that a state with C = is separable with p(q, q') = Pi(<z)p2 (<?')> an d 
thus a pure state. Indeed, this is a necessary criterion, and for the Weyl formulation this 
is equivalent to det T = 1. 
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4.3 The Klein-Gordon field 



As our main example of a bosonic field we use the simplest possible, namely the free Klein- 
Gordon field placed on a 1 + 1 dimensional chain with periodic boundary conditions. That 
is, we use the Lagrangian 



summed over all lattice points. The Euler-Lagrange equations gives the equations of motion 
and the dispersion relation 

u l = — sm 2 {k/2) + k 2 (4.18) 

where a is the lattice constant. To maintain a conformal invariance of the model, one must 
keep Na constant when rescaling the theory. Now follow the expectation value matrices 
for the ground state of the field, 

Q™ = ^E^' t(m - n} («») 



k 



P, 



mn ^ 



l -^k^ k(m ^ n) (4.20) 



k 



Smn = 0. (4.21) 



The sums over k are over all wave numbers 2nn/N with n = {0, ■ • ■ , N — 1}, discretized 
through the boundary conditions. 

In the massless limit, the field is expected to be conformal, and thus follows the conformal 
signature derived by Holzhey, Larsen, and Wilczek [HLW94] ISkr05aj (see Section 15.21) . 
Massless here means roughly k < 0.1. The true massless system k = is not accessible 
as Q diverges due to the k = term. A massive theory in the limit k — > oo carries no 
entanglement. 

One can also impose anti-periodic boundary conditions on the field, thereby explicitly vio- 
lating the translational invariance in the target space, and hence the conformal symmetry. 
The entanglement diverges as iV — > oo, like in the periodic case. However, in the limit 
k — > the periodic wave function leads to diverging entanglement, while the anti-periodic 
one does not. This is an example of entanglement divergence even though we do not have 
conformal invariance. Usually, in non-critical systems, the entanglement saturates at some 
point. 
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4.4 Two-mode entanglement 



The reduced density matrix of two sites in the lattice is generally a mixed state, and 
generally the entanglement in this state is ill-defined. However, for Gaussian CV states, it 
is known how to find the entanglement of formation for this state [ WGK + 04] . That is, the 
so-called Gaussian entanglement of formation Eq of a state pr with correlation matrix F 
is 



E (pr) = m£{S(pr,)\r p <r} 

1 v 



(4.22) 



where T p are correlation matrices of all pure states and S(pr p ) is the von Neumann entropy 
of the state. For the case of Gaussian states this is determined by the correlation matrix' 
symplectic eigenvalues. Given a symmetric, positive matrix T, one can find a symplectic 
matrix M G Sp(2iV) and a diagonal matrix D G R, Ar><Ar such that 

MTM T = D®D. (4.23) 

This is known as the symplectic diagonalization of T, and it is unique up to permutations. 
The elements of D = diag(di, • • • , d^) are known as the symplectic eigenvalues of T. Now, 
the von Neumann entropy of the pure state T p is |WGK+04j 



S(pi 



£c(«*) 



(4.24) 



where 



C(x) = cosh 2 (x) log 2 (cosh 2 x) — sinh 2 (x) log 2 (sinh 2 x) 
dk = coshafc, a>k > 0. 



(4.25) 



Given a 1+1D chain of bosonic systems, one can write the correlation matrix for two sites, 
known as the two-mode correlation matrix, as 



go qi 

p px 

<Jx Qo 

pi po 



"1 


0" 




"0 


0" 











1 



(4.26) 



where 

T q = [ go qi ] and T p = h Pl 

Here go = Qoo, Qi = Qon, Po = -Poo, and pi = Pan with n the distance between the two sites 
under consideration. Generally, any two-mode correlation matrix can under local unitary 
transformations be written in the so-called standard form 



n a 
k q 
k n 



k n 











n b 
n b 



(4.27) 
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with k q > \k p \. For the special case that n a = = n, the state is known as a symmetric 
state [GWK+03] and a necessary condition for this to be a correlation matrix is that 



n 



— kq > 1/4. Also, the state is entangled if and only if (n — k q ){n + k p ) < 1 [GKLCOl] 
Symmetric states arise for example as two light beams from a parametric down converter 
is sent through optical fibres or when a squeezed state is sent through two identical lossy 
fibres. The entanglement of the symmetric Gaussian state is particularly simple, 

E G (pr) = C(P), P = \ In [(n - k q )(n + k p )} (4.28) 

with ((x) defined in (14.251) . In terms of the correlation matrix this can be written as a 
symmetric state through the symplectic transformation 

S = diag (07, 1/ 07, 07, 1/ 07) 



where r\ is the squeezing factor 77 = ^/po/qo- This is a local, unitary Gaussian operation. 
Hence 



™ = 0oR) k q = q H — k p =p H —, (4.29) 

V Qo V Po 

and the entanglement follows from this. 

Some results for the Klein-Gordon field is shown in Fig. 14. 1L We see that the entanglement 
in the massless field vanishes quickly with distance, indeed for large systems the entangle- 
ment extends only to the nearest neighbor. For small systems it may extend to the third 
neighbor, due to the periodic boundary conditions. Also, all entanglement vanish for a 
massive system of k > 1, as is the case for the entanglement entropy [Skr05a] . However, 
there is no divergent entanglement for the two-mode case as opposed to the entanglement 
entropy, but the entropy saturates for iV — > 00, k — > at E G w 0.48. 



4.5 State evolution 

The massless Klein-Gordon field, or the bosonic vacuum investigated in [Skr05a] gives 
directly rise to a conformal symmetry and the properties that can be deduced from that. 
However, the time evolution of a state in vacuum is trivial, and to provoke a non-trivial 
time evolution of the system, we introduce an impurity and consider the toy model with 
Lagrangian density 

£n = I K - ( V ^™) 2 - «Vn] + ^nl- (4.30) 

This is the usual Klein-Gordon field bracketed and an impurity of strength e at a site / 
to provoke a non-trivial time evolution of the state. The system is still linear, and the 
discussion of this system is an example of how linear systems can still provide interesting 
behaviour. We assume the initial condition that the initial state is that of the unperturbed 
case e = 0. 
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Figure 4.1: The two- mode Gaussian entanglement between two sites in the ID Klein-Gordon 
field with periodic boundary conditions. 

Left: Eq for the nearest neighbor and next-nearest neighbor as function of total system size 
N. The data are for a massless (k = 10 -3 , x) and a massive (k = 10~°- 2 , □) system. The 
next-neighbor entanglement vanishes for large systems, while the 3rd neighbor entanglement is 
non-zero only for N = 6. 

Right: Eq for the three nearest neighbor pairs as function of n for a large (N = 30, x) and a 
small (N = 6, □) system. 



We apply periodic boundary conditions, and Fourier expanding the field in the Heisenberg 
picture gives 

with the initial condition eifc(0) = a k - The time evolution of the operators a and can 
now be computed by the Euler-Lagrange equations to be 

a k (t) = a k e^ + % (l - e"^*) , (4.32) 



where 



4 • 2, , /n\ i 2 , , I U k _ i/,/ 



uj k = — sin (k/2) + k and h k = \ ——se 
a 1 V 2iV 



Now, we consider the vacuum state that is the ground state of the unperturbed system 
where a k \^l) = 0. Then the correlation matrix is unchanged, but the expectation values of 
the field become 

(Vn(*)) = E ~2 t cos k ( n ~ l ) ~ cos - k(n - I))] , (4.33) 
and, if we consider the site of the impurity, n = I we find the equation of motion 

(<Pn(t)) = ^ E \ (! " 008 • ( 4 - 34 ) 
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When iV — ► oo we can write the constant term as 

l v - - £ r 1 

N 2-^< ,,,2 _ tt / ,,,2 ~~ 



7ra 



K 



V4 + 



2^2 



(4.35) 



Figures I4~2l show the time evolution of the field expectation value ((p n (t)) and its derivative, 
the expectation value of the conjugate field (it n (t)) = (<p n (t)). The field has a roughly 
cosine shape, governed by the main contributing k = term ("the zero-mode") of the sum 
14.331 Hence the period is 2tik~ 1 . Also, given that we consider the impurity term, a new 
perturbation hits this site every time a perturbation has propagated through the system 
to the site. Thus, with units where the perturbation speed in the system equals the lattice 
constant per time unit, there will be a new perturbation every N time units. If we consider 
another site 8 sites away from the perturbation site, the perturbation will hit every 5 + Nn 
and N — 5 + Nn time units with integer n. In the following we consider the impurity points 
and units where the perturbation speed is one lattice point per time unit and an invariant 
lattice constant a — 1. Now, we have two different regimes dependent on kN ^ 27r. For 
large systems where kN > 2tt, the main contribution is the oscillations of the point itself. 
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Figure 4.2: The field expectation values {<fi(t)) and (717 (t)) at the perturbation site I with time. 
To the left k = 1CT 2 and N = 1000 (kN = 10 > 2vr), and to the right k = 10 -3 and N = 500 
(kN = 0.5 < 2ir). Note that a perturbation hits every N time units and that the predominant 
oscillations have period 2itk~ 1 . 
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The field tends to relax its oscillation over time until a perturbation hits which excites 
the field again, followed by another relaxation period. For small systems with kN < 2tt 
the relaxation never occurs since the perturbation hits the point too quickly, and thus the 
perturbations are the main contribution to the time evolution. Figures 14.21 illustrate these 
two regimes. In both cases the leading oscillatory term is with frequency 2ttk~ 1 , even 
though the main oscillations have different origins. Also, the expectation value (iri(t)) has 
small oscillations of frequency 2tt, which in the figure appears as a smearing of the graph. 



5 Quantum Critical Systems 



Classical phase transitions occur under the change of temperature where at some critical 
temperature T c the system has an abrupt thermodynamical change, the hallmark example 
being the liquid-solid transition in e.g. water. Quantum phase transitions |Sac99j occur at 
zero temperature, and thus quantum effects can automatically be assumed to be profound. 
Indeed, the thermal fluctuation that take a classical state across a phase transition are 
non-existent, and the only fluctuations present are the quantum fluctuations determined 
by Heisenberg's uncertainty principle. 

The general setup is to consider a quantum system under the change of some external 
parameter 7, which causes a transition at some critical value 7 C , where a previously excited 
state now becomes the ground state. This crossover is called the quantum critical point. 
At this critical point there are, just as in classical phase transitions, scaling laws. For 
example the energy gap A, which necessarily is zero at the transition itself, will scale as 

A~ J\i-ic\ zv 

where zv is known as the critical exponent of the phase transition and J is some energy 
scale. Moreover, the correlation length £ will diverge (in an infinite system) at the critical 
point, and scale according to £ ~ L\j — 7 C | !/ . This also means that we have a scaling law 
between the correlation length and the energy gap, 

A ~ C z - 

Since the correlation length diverges - or extends to the entire system in the case of a finite 
system size - the system is scale invariant at the critical point, and the properties of the 
system becomes independent of the scale on which we view the system. 

Any classical phase transition invariably has temperature fluctuations which drive the phase 
transition. However, as the quantum phase transitions occur at T = 0, there can be no 
classical fluctuations in these models. Nevertheless there are always quantum fluctuations 
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Figure 5.1: Two scenarios for a quantum phase transition in which a parameter 7 is taken 
across the critical value 7 C . To the left there is no classical phase transition close to the quantum 
transition, while to the right there is a classical phase transition across the line, and in the shaded 
region the theory of classical phase transitions can be applied. [Sac99j 

due to Heisenberg's uncertainty principle, and these fluctuations have to drive the quantum 
transition, in the sense that they are necessary to explore the low-lying excitations of the 
system. Moreover, the classical correlations present in the classical models can be partially 
replaced by the quantum correlations which present themselves as entanglement. 

The quantum phase transition may or may not have a classical counterpart in the T > 
regime, see Fig. I5.1| in which case there will be a region where there are competing classical 
and quantum fluctuations across the boundary. However, there are also cases where the 
thermodynamic variables are analytical close to the T = line, and thus no classical phase 
transition takes place. 

In the study of quantum statistical mechanics, we use the density operator, which at finite 
temperatures take the form 

P = e-** 

where (3 is the inverse temperature and H the Hamiltonian. However, in the zero-temperature 
limit (3 — > 00 only the ground state will contribute to this expression, and assuming the 
Hamiltonian is diagonal in some basis H = Y2 n E n \n){n\ with lowest energy level Eq = 0, 
the density operator becomes 

p (T = 0) = e-^En\n)(n\ = e -/3£o|0><0| = [ > < [ . 

Hence, when computing properties of a quantum critical system, we need to identify the 
ground state and hence follows the density matrix and all thermodynamic properties. 

A phase transition can, strictly speaking, only occur in the thermodynamic limit of ~ 10 23 
particles, which is naturally inaccessible when implementing the system on a computer. 
However, one may see traces of what may become a phase transition at smaller systems, 
even with as little as ~ 10 particles [SO05] . However, it is important to emphasize that no 
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actual phase transition can occur in this system, and that when we speak about parameters 
where we get a phase transition, it should be interpreted as parameters where a phase 
transition occurs in the thermodynamic limit. 



5.1 Ground state entropy 

When the ground state of a quantum system changes macroscopic features across a quan- 
tum phase transition, the details of the ground state is vastly more complicated at the 
phase transition point than at non-critical places in the parameter space. This is indicated 
by the increase of entanglement at the critical point as measured both by the concurrence 
[OAFF02] and the entanglement entropy [VLRK021 ILRV04] . As the entropy involves trac- 
ing out some part of the system, it will have a dependence on the size of the system traced 
out, as well as the overall system size. However, it should not depend on the small-scale de- 
tails of the models since these are irrelevant at critical points according to renormalization 
theory. 

Specifically, at the critical point in theories where one has conformal invariance in 1 + 1 
dimension, the entropy diverges logarithmically with system size, as opposed to off-critical 
points where the entropy saturates [LRV04j . However, close to the critical point, the system 
size where the entropy saturates is very large, and thus this is not a very efficient indicator 
of whether or not a system is critical due to the computational resources required. 



5.2 Consequences of conformal symmetry 

omce Polyakov |Pol70] it has been known that local scale invariance implies conformal 
invariance and that this fact can be used to investigate criticality in statistical mechanical 
systems. Belavin, Ployakov and Zamolodchikov [BPZ84J showed how the rich structure 
of two-dimensional conformal theories can be used to extract information on the critical 
properties of a quantum field theory. 

A quantum field theory becomes conformal if it fulfills the following requirements [FMS97, 
IGin89j : 

• Translation and rotation invariance 

• Local interactions 

• Lorentz invariance 

• Scale invariance 
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The last condition is the crucial one, since it is only fulfilled at critical points. Thus we 
expect conformal symmetry on the critical point, but not elsewhere in the parameter space. 
The conformal group in two dimensions is infinite as opposed to larger dimension^. That is, 
any holomorphic transformation is conformal, and to specify a holomorphic transformation 
means to know all the coefficients of its Laurent series, which are infinitely many. This 
means that conformal symmetries in two dimensions are easier to analyze and has a richer 
structure than those in larger dimensionalities, and this is the reason that we will exclusively 
focus on this case. 

We define the real coordinates (z°, z 1 ) on the plane, and introduce complex coordinates, 

z = z° + iz 1 z = z° — iz , 

with differentiation rules 

d = d z = -(<9 -idi) 9 = dg = -(<9 + i<9i). 
Under a holomorphic transformation z i— > f(z), a line element ds 2 = dzdz transforms as 

— (!)(!>> 

This transformation law is generalized under the assumption that $(2, z) dz h dz h is invari- 
ant under the transformation. Here Q(z,z) is a field and h and h are real parameters. 
Then the field transforms as 

$( ^^ (f )\§y $ ( /(z)j( ^- (s - i} 

This defines a primary field of conformal weight (h, h). Non-primary fields are known as 
secondary fields. 

A key ingredient when computing correlators in conformal symmetry is the operator prod- 
uct expansion (OPE) . This is the assumption that the product of two local fields 4>i{ x ) 
and 4>j{y) can be written as a linear combination of local operators, 

<t>i{x)<t>Ay) = Y, C U X - v)Mv) 

k 

where the coefficients are c-numbers. Now, one is usually interested in the singular 
behavior as x — > y, and we use ~ to indicate that regular terms are discarded. For a 
primary field (f with conformal dimension h, the OPE of this with the energy momentum 
tensor is ^ ^ 

T(z)(p(w,w) ~ -=Lp(w,w) H d w Lp(w,w) 

(z — wr z — w 



The one dimensional case is trivial, since any continuous transformation is conformal. 
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and correspondingly for the anti-holomorphic fields. For the energy-momentum tensor with 
itself the OPE is 

Tl W \ 0/2 2T( W ) &T(w) 

T{z)T{w) ~ - + - - + . 5.2 

(z — wp [z — w) 2 z — w 

The constant c is called the central charge of the specific model in question, or sometimes 
denoted the conformal anomaly. Apart from this anomaly, the energy-momentum tensor 
is simply a conformal field of conformal dimension h = 2. Of course, the relation (15.2ft 
holds also for the anti-holomorphic parts of the theory, thereby defining an independent 
anti-holomorphic central charge c. When the theory has a Lorentz-invariant and conserved 
two-point function (T^ u (z)T a p(—z)) , one must have c = c [Gin89|. This holds for all 
relevant quantum statistical systems. 

Central in conformal field theories is the study of correlation functions, a subject with 
obvious connections to critical systems. In particular, the two-point correlation function 
of two fields 0i and <p 2 is nonzero if and only if they have the same conformal weights 
hi = h 2 — h and hi — h 2 = h. Then, 

(0l(Zi,Zi)0 2 (z 2 ,22)) 



(zi - z 2 y h ( Zl - z 2 f h 

Notably, the correlation function depends only on the distances Z\ — z 2 and Zi — z 2 , a 
statement that is true also for the three- and four-point correlation functions. 

The Laurent expansion of the energy-momentum tensor can be written 

oo oo 

T(z) = z ~ Z ' 2L n and f{z) = £ z~ z ~ 2 L n , (5.3) 



n=— oo 



and the inverse relations are 



dZ z n+1 T(z) and L n = I ^ z n+1 f (z). (5.4) 



2vri v 7 " / 27ri 

Now, one can compute the commutators of L n by performing the integrals and expressing 
the OPE of the energy-momentum tensor, and the results are 

Q 3 

[Lni L m ] — (jl — m)L n+m + ~^{ n ~ n )$n-m 

— — — c 

[L n , L m ] = (n - m)L n+m + — (n 3 - n)5 n - m 
[L n , L m ] = 0. 

Thus we have two independent algebras on L n and L n of infinite dimension. This is known 
as the Virasoro algebra of the central charge c. If c = 0, the algebra defined on L n is the 
classical algebra of the generators for infinitesimal conformal transformations in the plane. 
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5.2.1 The free boson 



The free massless boson in two dimensions has an Euclidean action [FMS97J 

S[<p] = ±g J d 2 xd&d»<p, (5.5) 

where g is a normalization constant. The two-point correlator of this field is 

((p(x)<f(y)) = -— ln(x - yf + const. 

When introducing the complex coordinates z and z, this becomes 

(<p(z, z)ip(w, w)) = [ln(z — w) + ln(z — w)] + const, 

4vrc/ 

or when separating the holomorphic and anti-holomorphic parts through a differentiation, 
the operator product expansion of the holomorphic field dip = d z ip becomes 

(d z <p(z,z)d w ip(w,w)) = 



Aug (z — w) 2 
1 1 

d(f(z)d(p(w) ~ -- — ^. 

4irg [z — w) z 

The bosonic symmetry property of this OPE is immediately obvious. Now, the energy 
momentum tensor in complex coordinates is 

T(z) = -2irg : dip(z)dip{z) : (5.6) 

Here : • : denotes normal ordering to ensure that vacuum expectation values vanish. Wick's 
theorem allows us to compute the OPE of the energy-momentum tensor with itself; 

T(z)T(w) = 47rV : dtp(z)d<p(z) :: d<p(w)d<p(w) : 



37r 2 g 2 : d(f(z)dip(z) :: dtp(w)dtp(w) : 



167T g : dif(z)d(p(z) :: d<p(w)dip(w) : 

1/2 : dif{z)dip{w) : 

Aixg- 



(z — wY " (z — w) 2 
1/2 2T(w) dT(w) 



(z — w) 4 (z — w) 2 (z — w) 

the last equation arises by a linear expansion around w, dip(z) = d<p(w) + (z — w)d 2 ip(w). 
Thus, the OPE has an anomalous term (1/2) /(z — w) 4 which fixes the central charge of 
the model at c = 1, and the energy- momentum tensor is not a primary field. 
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5.2.2 The free fermion 

The free fermion in two dimensions is described by the two-component spinor \1/ = (if),ip), 
and the action is 

S(V)=g [ d 2 x(tpdtp + tpdtp). (5.7) 



The propagator of this field becomes 

1 1 



(il>(z,z)ip(w,w)) 



2irg z — w 
and thus the OPE follows trivially 

i/)(z)ip(w) 



2ng z — w 

which of course has fermionic symmetry properties under particle exchange. The energy- 
momentum tensor is now 

T{z) = -ng : ^(z)d^(z) :, 
and by Wick's theorem we arrive at the OPE for T(z) with itself; 

T(z)T(w) = 7rV : ip(z)dip(z) :: ip(w)d^(w) : 

1/4 ng 1 

+ — - rr (: ip(z)ip(w) : — : ip(z)dip(w) : +dif)(z)if)(w) :) 



[z-wY 2 (z-w) 2 
~ — " — : dtp (z) dip (w) : 
1/4 2T(w) dT(w) 



+ 7 77T + 



[z — w) 4 (z — w) 2 z — w 



thus proving that the free fermion has a central charge c = 1/2. 



5.2.3 Entanglement entropy in a 1 + 1 dimensional strip 

After Hawking's discovery of black hole radiation in 1975 |Haw75] . the study of quantum 
effects and information exchange in black holes became a topic of intense study, and thus 
followed investigations on entropy in conformal field theory that have resurfaced for use 
in quantum information theory. Srednicki |5re93j computed numerically the entropy of 
both 3 + 1 and 1 + 1 dimensional bosonic, massless models and identified the logarithmic 
divergence of the latter. Later, Callan and Wilzcek |CW94] came up with the concept of 
geometric entropy on a 1 + 1 dimensional conformal strip, which was pursued by Holzhey, 
Larsen and Wilczek [HLW94J. The results have been generalized by Calabrese and Cardy 
to other geometries and finite temperature |CC04l IKor04j . and also the time evolution of 
the entanglement has been investigated [CC05] . 
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The analysis starts with a series of conformal transformations in the plane. Assume that 
the initial coordinate is z = x + ir, where x is the spatial coordinate and r is the time. 
The system is of length L, such that < x < L, and the subsystem we're interested in 
is < x < i. We denote the surfaces at r = as S for the subsystem and C for the 
entire system. Now, periodic boundary condition in x make this model an infinitely long 
cylinder of circumference L, the length of the system. The first conformal transformation 
is |HLW94j 

sin If (z -£)] 

C = L . £ Lx • (5-8) 

sm (t) 

This maps C onto the real axis, with S as the negative half- axis and the rest S' = C\S as the 
positive axis. Moreover, the infinite past is a point in the upper half-plane, £(r = — oo) = 
_ e -mi/L^ anc j i n fi n ite future is at the point £(r = oo) = — e mt / L in the lower half-plane. 
Next, we employ the transformation w = ~ln£, which maps S onto the real axis with 
x — > corresponding to w — > oo, and x — > £~ corresponding to w ^ — oo. The remainder 
of the system is the line w = i with z — > £ + becoming w — > — oo + i, and z — > L becoming 
w — >• 00 + i. The infinite past and future maps to the points w{r = ±00) = i(l ±£/L). 
Both transformations z — > C — > w are sketched in Figure | 



Now, we can use the so-called replica trick to make the identification 

S = — Tr p log p = — lim -^-Tr p n . (5.9) 
rwi cm 

This assumes that integer values of n can be made into an analytical continuation for 
non-integer n. Thus the problem is reduced into computing p n on the surface defined by 
w. The functional integral defining p is 

p(0,0') = J 

where \ designates the fields on the upper line, while <fi are the fields on S, that is the real 
axis. This integral is equivalent to cloning the strips C, with the intersection along S', with 
boundary conditions specifying <p and <$' . Thus one gets a strip of width 2 to integrate 
along. Using functional integral representation of the wave function \<j),x) one finds that 



where Z{n) is the path integral of the strip C cloned n times such that one gets a structure 
of width 2n and boundary conditions restricting the fields on the extremal surfaces. From 
this one finds the entropy (or entanglement entropy) to be [HLW94, CC04J 



c + c 
o ~ — - — log 



(5.10) 
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Figure 5.2: The conformal transformations described in the text, z 
( —* w = ^lnC- The red arc indicates the system, S. 



C defined by (|5.8p and 



where ~ indicates up to an additive constant independent of i and L. The derivation is 
valid in the limit £,L — i 3> a with a the lattice spacing of the model. The entropy is 
symmetric around I = L/2, at which point the entropy is largest. From this it is easy 
to see that for an (semi-)infinite system, L — ► oo, the entropy diverges logarithmically 
S ~ g -rr L \ogi. Also, when considering systems of different sizes L, with a constant fraction 



/ L, one also has a logarithmic divergence, S 



log Zr. 



5.2.4 Conformal invariance in QPT 



We expect conformal invariance at quantum critical points in 1 + 1 dimensional models due 
to the scale invariance that emerges here. It is reasonable to demand nearest-neighbor (or 
perhaps next nearest-neighbor) interaction to keep interactions local. Thus it is possible to 
identify a central charge to the quantum phase transition (QPT), and the central charge 
defines the corresponding Virasoro algebra. Hence it is, in principle, possible to identify the 
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critical exponents of the model according to the scaling dimensions of the fields. Friedan et 
al. |FQS84| proved that in models where the Virasoro algebra has unitary representation^], 
and c < 1, the central charge is restricted to the values 

6 

c=l-— r — - m = 3,4,---, (5.11) 
mym + 1) 

and for each value of c there are mim — l)/2 allowed values of h given by 

[(m + l)r — ms] 2 — 1 
h rs (m) = — — r = 1, • • • , m — 1 and s = 1, •••,50. 

' K ' 4m(m + l) ' 

Moreover, each such value of c corresponds to a universality class (or a set of such), and 
for the first few values the mapping is shown in Table 15.11 

Table 5.1: The universality classes corresponding to the first allowed central charges c < 1. 
Fqg84llUrn89] 



m c Universality class 

3 | Ising, free fermion 

4 jq Tricritical Ising 

5 | 3-state Potts 

6 | Tricritical 3-state Potts 



oo 1 Free boson 



It is remarkable that the conformal invariance at the critical point gives information about 
which universality class the model belongs to. There is a deep connection between the 
conformal invariance of the underlying structure and the statistical mechanics of the model, 
a feature that we will exploit in later chapters. 



5.3 Quantum models 



We will in this section present a limited number of quantum models which exhibits quantum 
phase transitions and describe some properties. In general, a spin-1/2 quantum chain with 
periodic boundary conditions has a Hamiltonian 



n=l 



fa &n a n+l + 9 ' (°"n X (T n+1 ) + A ■ (T n 



_a=x,y,z 



(5.12) 



2 Such models are usually called unitary models, and are those where the representation of the Virasoro 
algebra contains no negative norm states. All such representations with c > 1 (and h > 0) are unitary. 
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where /, g, and A are real 3-vectors, and <tn+i = &i- Further, cr n = (a^,a^,a^) are 
the conventional Pauli matrices. Fixing the coordinate system, this leaves us with seven 
degrees of freedom in the Hamiltonian, which is a huge parameter space to explore for 
phase transitions. The models with various parameters zero, or equal are usually known 
under a plethora of names. We will always consider g = 0, and assume A to point in the 
positive ^-direction, A = (0, 0, A). Still the model can be adjusted by a constant, so we can 
write down a general three parameter model, which we will denote the XYZ model; 



Hxyz — — 



with f x + f y = l. (5.13) 



When f y = f z = this is known as the quantum Ising model, and when f y > we get 
the XY model, both of which will be discussed in detail shortly. In particular the critical 
XY model with f x = f y and f z = A = is sometimes known as the XX model [WanOlb] . 
When f x = f y = f z > and A > we have the XXX model or simply the ID Heisenberg 
model (ABVOlj . When A = and / > we get what is also sometimes known as the 
XYZ model, a model that in terms of transfer matrices can be considered equivalent to 
the classical eight- vertex model |Sut70j. In the same way, the case A = 0, f x = f y > 
an d f z > is known as the XXZ model, or in terms of the transfer matrix equivalent 
to the classical six- vertex model [Lie67l IWanOla] . which in its quantum version has been 
suggested to describe the <i-density wave competing with superconductivity in high-T c 
superconductors [Cha02t [SC06J. Curiously, the model was in its early days considered 



for explaining "residual entropy" in ice at low temperatures |Lie67j . which suggests the 
wide applicability of these models, not restricted to quantum mechanics. All these models 
investigated and some solved exactly, which makes them useful as benchmark tools. In 
particular, the Ising model is simple and intuitively easy to assess. The various models are 
not only mathematical playgrounds - though they are useful as such - but are also efficient 
models to simulate real life systems in addition to the examples already mentioned. E.g. 
the vertex models were originally conceived to model a second order polarization transition 
at 122 K in potassium dihydrogen phosphate, KH2PO4 |Sla41j . 



5.3.1 The Ising model 

The Ising model is perhaps the quantum model which is studied in greatest detail, and 
where most is known. Hence it makes for a good starting point for further studies. We 
investigate the quantum Ising model in one spatial dimension, whose Hamiltonian is 

H Ising = -X>x+i+A<). (5 - 14) 

n=l 

The as are the conventional Pauli spin matrices while A is our external parameter, a 
magnetic field in the z-direction. Naively, one expects in the limit A — > 00 all spins aligned 
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in the z-direction, which is a product state and thus has no entanglement. However, in 
the zero field limit, the ground state would have to be a Schrodinger cat state with a 
superposition between the spins in positive and negative x-direction respectively. 

This ID quantum model can be mapped onto the 2D classical Ising model, which Onsager 
famously solved in 1944 [Ons44] . The classical model has a continuous phase transition 
at a critical temperature T = T c , which in the quantum case is mapped onto a critical 
magnetic field A c = 1. 

In the A — > limit we have unit entropy due to the one-qubit entanglement in the ground 
state, which is a Schrodinger cat state, while the large A limit the product state means 
that the entropy falls to zero. However, we see that with larger system sizes there appears 
what seems to be a divergence at A = 1. Even in small systems the entropy here is larger 
than elsewhere, but for larger systems the effect is more pronounced. Note that we need a 
large system (roughly iV > 50) to say with any confidence that there is a phase transition 
at A = 1, which indeed is the thermodynamic solution to this system. 

In Figure l5\3l the gap and ground state energy of the Ising model is shown, and the approach 
towards a gapless state as iV — > oo is clearly shown. Note also that the N = 10 system is far 
from gapless, and one should therefore not anticipate much sign of the phase transition in 
this system, contrary to what shown in [SO05] and Section 15.4.51 Indeed, it is remarkable 
that even though a system of size of order 10 is far from gapless, and shows no critical 
properties, the conformal signature can still be very reliably used to identify the critical 




Figure 5.3: The gap AE = min n u) n between the ground state and first excited state energies 
in the Ising model with periodic boundary conditions. From top, the system sizes shown are 
N = 10, 50, 100. We see clearly that the system approaches a gapless state at the critical point 
A = 1 as the system size increases. The inset shows the ground state energy for the N = 100 
system. 
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parameters. The gap is here defined as AE = min n u; n , within the ground state parity 
segment. 

5.3.2 The XX model 

The XX model in an external magnetic field is defined by the Hamiltonian 

N 

H xx = -J2 («+i + «+i + A <) • ( 5 - 15 ) 

n=l 

This model has a continuous U(l) symmetry in the xy plane, which by the Mermin- 
Wagner theorem |MW66] cannot be spontaneously broken at non-zero temperatures in 
two dimensions in classical systems. Thus, in quantum models in one dimension at zero 
temperature, no spontaneous symmetry breaking can occur. Nevertheless, polynomial 
decay of the correlation function may exist due to vortex excitations such as the Kosterlitz- 
Thouless transition in the classical XY model |KT73j . The parity operator V = <S>k°~k 
commutes with the Hamiltonian, and thus its eigenvalue, the parity V = ±1 is a good 
quantum number, and the solution space is split into two subspaces, as demonstrated in 
Section 13.11 In this model the generator of rotation, C = J2 n a n a l so commutes with the 
Hamiltonian, and its eigenvalue I is thus also a good quantum number. In terms of the 
spin flip operators a ± = (a x ± ia y ) / \/2, the Hamiltonian can be rewritten 

H X X = ~J2 Vn+l + CT n ff n+l) ~ X ^ 

n 

Obviously, the product states \N) = \l) m and | - N) = \0) m (labelled by their / values) 
are the exact ground states in the limits A — > ±oo respectively. More precisely, applying 
the Hamiltonian to these states yield H X x\N) = —XN\N) and H X x\ - N) = XN\N), 
showing that the state \N) has lower energy of these two when A > 0. However, it is not 
obvious that these states are the ground states for all A. The parities of the states are 
p(\N)) = 1 and p(\ - N)) = {-1) N . 

Consider the simplest "excitations" to the state \N), the iV states 

\N - 2, m) = 1 1)®^- 1 ) ® |0) ® 1 1)®^-" 1 ) m = 1, • • • ,N. 

Applying the Hamiltonian to this state gives 

H xx \2-N,m) = -[\N- 2,m- 1) + |iV- 2,m + 1) + X(N - 2)|2 - N,m)} . 

Thus, it is not an eigenstate of the Hamiltonian, but we can define a translationally invari- 
ant state, or a spin-wave state 

|iV-2,B = \"e 2wikm/N \N -2,m), fc = 1, ■ ■ • ,N, 
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which diagonalizes the Hamiltonian. That is, 

N-2,k}. (5.16) 

The excitation energy of the spin- wave state compared to the product state \N) is AE = 
+2 (A — cos ^jf), which may be negative when A < 1. Likewise, we can define the simplest 
excitation to the state | — N) as 

| -N + 2,m) = \0)®( m -V g, |i) g, \ty®(N- m ) m — 1, ■ ■ ■ , N, 

and the corresponding spin- wave state | — N + 2, k}, which diagonalizes the Hamiltonian 
and has energy eigenvalues 

-N + 2,k}. (5.17) 

Here, the excitation energy is AE = —2 (A + cos ^fr), which can be negative when A > — 1. 
Thus, the ground state is shown to be \N) when A > 1 and | — N) when A < — 1, while the 
ground state in the intermediate region is undetermined. 



H xx \N-2,k} 



-XN + 2 A 



cos • 



2nk 
~N~ 



H xx \-N + 2,k} 



XN — 2 A + cos 



2nk 



5.4 The XY model 

The XY model is a simple generalization of the Ising and XX models, conventionally 
written 

Hxy = -J2 (^r^X+i + + A<) , (5-18) 

where 7 = 1 is the Ising model and 7 = the XX model already discussed. First 
properly solved by Barouch and McCoy [BM7TJ this is now a thoroughly investigated 
model |Ort05t ION02a] which is useful for benchmarking due to its simplicity and still rich 
structure. A schematic phase diagram of the model is presented in Figure 15.41 A simple 
assessment of the model shows that the ground state in the large field limit where A — ► 00 
is the pure state \Q Z ) = \ |) 0Ar . In the low field, when 7 > 0, the ground state will be 
the Schrodinger cat state with positive and negative x direction being equivalent, \Q X ) = 
(I -^.}® N -f J /v2. Finally, for negative 7 the ground state will be the Schrodinger 

cat state with superpositions in the y direction, \Q y ) = (| y}® N + | y)® N ^} j\f2. Here 
vertical arrows indicate eigenstates of a z , horizontal arrows eigenstates of a x and diagonal 
arrows eigenstates of a y . Strictly speaking these assessments are only valid in the extreme 
limits, but one can assume that they are reasonably close to the actual states also in some 
region away from these limits. In the thermodynamic limit one can assume this to be the 
case arbitrarily close to the phase transition. 
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Figure 5.4: The phase diagram of the XY model as defined in (|5.18j) . The blue lines indicates the 
phase transition between the pure state with spins aligned in the z direction and the Schrodinger 
cat states. The red separates the Schrodinger cat state in ±x direction from ±y direction. The 
circle indicates the BM-circle. 

Barouch and McCoy also acknowledged the existence of the unit circle 7 2 + A 2 = 1 as 
the division line between what they denoted the oscillatory region inside the circle and 
monotonic region outside due to the behavior of the xx correlation function in the large 
distance limit |BM71j . However, this is not considered a true phase transition, but rather 
a boundary region. The three phases of the model have been denoted ordered oscillatory 
(7 2 + A 2 < 1), ordered ferromagnetic (7 2 + A 2 > 1 and |A| < 1) and paramagnetic (|A| > 1) 
[WDM + 05| . At this so-called BM-circle after its discoverers, the entanglement entropy is 
always unity [LLRV05J. 



5.4.1 Fermionization 

One of the reasons for the feasibility of the XY model is that it can be mapped onto a string 
of spinless fermions, using the technique of a Jordan- Wigner transform. This reduces the 
dimensionality of the problem from the 2 N x 2 N matrices of the spin formulation to 2N x 2N 
sized matrices. The technique can be applied more generally, but we will demonstrate it 
on the XY model here. 

We apply a Jordan- Wigner transform that introduces fermionic operators a n and 




(5.19) 
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which ensures the anti-commutation of these operators, 

The spin operators can be expressed in terms of these fermion operators as 

-t -1 



i 2(l n Q Jn 



«+l 




"-1 


-1 


1 1 " 


°>l+l 




i 


— i 


— i i 


«+l 




i 


i 


i i 


_°n°n+l_ 




1 


-1 


1 -1 



. a n a n+l. 



(5.20) 

(5.21) 



In turn, this means that the Hamiltonian can be written as a sum of quadratic fermion 
operators, which in general is 

with A and B being N x N matrices where A is real and symmetric. In the special case 
of the XY model these parameter matrices become 



.4 



A 

_i 

2 





A 



B 



2 7 



0-10 
1 
1 
1 



1 



(5.22) 



The above holds for the case of open boundary conditions, where the sum in (I5.18P goes 
from 1 to N — 1. However, for the case of periodic boundary conditions, the upper limit 
of the sum is N, and by definition <Jn+i = &i- For this case we must utilize the parity 
operator V = <^)^ =1 a^, which commutes with the Hamiltonian, and a^ai = —a^a^V. □ 
Also, for generality one can assume that the chain is non-homogeneous, such that the A 
and 7s are local, {A n } and {7 n }- For the case of the A and B matrices above, they now 
become 



.4 



Ai 




\v 







1 X 1 

2 A 2 -5 
"2 ^3 
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-7i 
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• - lN V 
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• 7iV-i 



3 To see this, use direct computation and the fact that tr d 



T(T ± CT Z . 



(5.23) 
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Replacing homogeneous parameters and a virtual parity V = 0, we recover Eq. (15.221) . 
Hence, the open boundary conditions can easily be recovered from the equations for periodic 
boundary conditions by setting V — 0. 

For the homogeneous, closed chain we thus obtain 

1 N 

H X y = - ^2 [ a " a «+i _ a l a n+i + 7 (a n a n+ i - alal +1 ^ + A (a) n a n - a n a) n ) (5.24) 

n=l 

where we must interpret clncln+i = — aN&\P ■ Now, considering the positive parity sub- 
space, V = 1, the boundary condition in terms of the fermions become anti-periodic, and 

we can expand the field in an anti periodic Fourier series of new fermionic operators b k and 
6 t. 



1 - 



iri(2k+l)n/N 
k" i 



1 - 



-wi(2k+l)n/N 



N 

r(+) 



k=l 



In terms of these, the positive parity Hamiltonian H X y becomes 

N r 



H 



(+) 

XY 



£ 

k=l 



COS 



7T(2k+ 1) 

N 



a (hbi - b\b k ) + i 7 



sin 



n(2k + 1) 
N 



b- k b k ~ b[bl k 



We have used the convention that —k = N — k — l. Next, we use a Bogoliubov transfor- 
mation to diagonalize the Hamiltonian, 



ib k 




-at*. 





cos <pk — sin <pk 
sin 0^ cos 4>k 



where the parameter is assumed to fulfill = —</>/. and Cfc,c£ are new fermion operators. 
For brevity, we define a k = ^2h±H ; anc [ g e t 

ffxy = { K cos ak ~ A) cos 2< ^ fc ~ 7 sin «fc sin 2 0fc] (c fc c|. - 4 c fc) 

+ [(cos ct fc - A) sin 2<p k + 7 sin a k cos 20 fc ] (c_ fc c fc + 4 c -&) } • 



This becomes diagonal upon choosing 

A — cos a k 



cos 20£ 



and 



sinz^fc 



7 s m «fc 

A4 ' 



with 



A/"fc = (A — cosafc) 2 + 7 2 sin 2 aifc. 

Note that the condition = — <fr k is fulfilled. This finally means that the Hamiltonian 

becomes 



N 

H xy = ^2\^~ cosa fc ) 2 + 7 2 sin 2 a k \ 2c\c k - l\ 

k=l s v ' 



(5.25) 
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The same procedure can be applied to the negative parity subspace, only now with periodic 
boundary conditions and periodic Fourier series, 



N 



2-rcikn/N 



k=l 



and we arrive at (I5.25P with (%%■ = InkjX and corresponding new fermion operators. In 
Figure 15.51 the energy eigenvalues is plotted with for various parameters, and we see 
clearly that the critical points correspond to zeros in the energy gap. 

A 




Qa 



Qi 



Q-2 



7 



Figure 5.5: Left: The energy modes £lk are shown for the four points in the phase space (7, A); 

Qa = (0, 0.5), Q 2 = (1, 0.5), Q 3 = (1, 1), and Q 4 = (0, 1). 

Right: The points shown designated in the phase space of the XY model. 



The parity of the ground state is important when investigating small systems — for large 
systems the difference is minute — in order to know which subspace to look for the ground 
state. This is not a trivial question, and we resort to direct computation of the two 
subspaces, and check which yields the lowest ground state energy. Brute computation 
shows that when 7 2 + A 2 > 1, that is outside the BM-circle, the ground state parity is 
always positivd3- Inside the BM-circle, the nature of the ground state is oscillatory, and 
the parity of the ground state in a N = 10 system is shown in Figure 15.61 The BM-circle 
limits the outer region of the circle, and the oscillatory nature of the ground state inside 
the circle is shown. When the system size increases the difference between the ground state 
in the different parity subspaces become negligible. 



5.4.2 Diagonalization in terms of Majorana fermions 

Having established the fermionic shape of the Hamiltonian in question, which in this case 
is the XY model, and which take on the form of the two matrices A and B, we must 

4 When N is odd, it will be negative outside this circle for negative A. We consider only even N in the 
following. 
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Figure 5.6: The areas in the first quadrant of the parameter space of the XY model where the 
ground state has negative parity are shaded for an N = 10 system. The outermost line coincides 
precisely with the BM circle 7 2 + A 2 = 1. The critical lines and the Ising line (dashed) is drawn 
for reference. 



diagonalize the Hamiltonian. To this end, define the 2N Majorana fermions in terms of 
the operators 



1 < t 
7w-i = ~/p a * ~ a i 



1 , u 
l2i = -j=(ai + al). 



(5.26) 



These operators satisfy the Majorana anti-commutation relations {7i,7j} = Moreover, 
in terms of these operators the Hamiltonian becomes 



where the Hermit ian matrix C is 



C = i 



-ii 

-21 



-12 
-22 



Im B mn 
A — Rp R 



' \ -^-mn ~t~ f^ 

Im B mn 



C being imaginary and antisymmetric means that its eigenvalues are real and come in 
positive/negative pairs, which we denote ±u) n . The eigenvectors corresponding to such a 
pair are complex conjugates. In turn we now diagonalize C 2 with the orthogonal and real 
matrix O, 



N 



o?c 2 o = 



n=l 



'"1 o 
ul 



50 



Entanglement and Its Applications in Systems with Many Degrees of Freedom 



If we are careful in permuting the columns of O, we can block diagonalize C such that 

N 

o?co = 



N 

—lUr. 



\uj n 



71=1 



Now we can return to fermions again, defining the N fermionic operators 



2N 



i=l 

These diagonalize the Hamiltonian, such that 

N 



a n = -i 7*(Oi,2n-i + iOi,2n) (5-27) 
^ 2 ,=i 

^ 2iV 

a « = 7/f S7i(Oi,2n-i - iO ii2n ). (5.28) 



if = w n (a+a„ - a„a^) . (5.29) 

n=l 

Hence it is clear that the eigenstates of the Hamiltonian will be states with a defined 
number of a-fermions. How this maps back to the real space fermions — never mind the 
spins — is far from trivial though. Nevertheless, it is clear that the ground state is the one 
with zero a fermions. Any state in this basis we denote (77) = |ni,ri2, • • • ,un) where r\ is 
the set of fermionic occupation numbers and — 0, 1. 

This shows how the orthogonal matrix O plays a crucial role in diagonalizing the Hamil- 
tonian, and that the new Hamiltonian is diagonal in terms of quasi-fermions related to the 
original fermions through the Bogoliubov transformation 

a m = S mn a„ + A m „4 a) m = A* mn a n + £ mn a m (5.30) 

where the transformation matrices are 

'-'ran n [0 2n ,2m- 1 + 2n -l,2m + 1 (0 2 n,2m ~ 2n -l,2m-l)] 
A mn = - [p2n,2m-l _ ^2n-l,2m + 1 {C>2n,2m + 02n-l,2m-l)] ■ 

Preservation of commutation relations and unitarity of the transformation ensure the fol- 
lowing properties of the matrices X = {S mri } and A = {A mn }, 

SSt + AA^l (5.31a) 

S t S + A T A* = l (5.31b) 

SA T + AS T = (5.31c) 

S t A + A T S* = (5.31d) 
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Hence, when transforming from real fermions to the quasi-fermions, we apply the columns 
of O, where for each particle labeled n, two columns of O contribute, and every second 
term in this column contributes at a time. There is as we see a complex structure of the 
orthogonal matrix that creates the quasi-fermions. In general, the ground state of the 
Hamiltonian (15.291) will be 1 000 • ■ ■ 0) in terms of the quasi-fermions, while it is not obvious 
how this can be mapped back to the real fermions. Nevertheless, this provides us with an 
easy manner in which we can find the properties of the ground state, such as the entropy. 

The following argument ensures that we are able to diagonalize the basis we are working 
in, based on a matrix notation of (15.30!) . 



a 




£ 


A 




a 






A* 


£* 




at 



Here, a is a column vector consisting of the a n 's, while a' is a column vector consisting of 
a^, not the adjoint of a. We can diagonalize the matrices EE* and AA^ simultaneously, 
since they commute by equation (I5.31al) . Hence, we can write 

^-UE%tf and AA f - rm2 



UD\U\ 



U being unitary and Dj^a are diagonal. Precisely likewise, we can do the same for the 
matrices involved in fl5.31bl) . 

£ t £ = V J D|V t and A T A* = (A* A)* = VD 2 A V\ 

V being another unitary matrix. The eigenvalues are however the same. Now, in order to 
preserve the two latter conditions in 15.311 we must have 

{D A ,D x } = 0. 

This can be accomplished since we have double degeneracy in both matrices by making 
these matrices 2x2 block diagonal with two different Pauli matrices along the diagonal, 

N/2 N/2 

^ E = COS ^ ff2 an d -Da 

k=l k=l 

These will be diagonal when squared, and thus comply with the conditions above. We 
have assumed that the two degenerate terms are coincidental, if not the matrices can be 
rearranged to conform with this. If N is odd, one will have to add a single term along 
the diagonal consisting of a (1) for the (-Da)- This procedure will ensure that the 
Bogoliubov transformation can be written 



0sin^ fc e iXk a y . 



->2n-l 
u 2n . 



COS$ n 

— i sin i? n e" lx " 



— i sin -# n e 1Xr ' 

— COS $ n 



i>2n-l 

bt, 



J 2n 



where the operators are unitarily transformed, 



kQ-k- 



(5.32) 



(5.33) 



Hence we have shown that the diagonalization can be accomplished through the unitary 
transformation (15.331) and Bogoliubov transformation (15.321) . 
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5.4.3 Entropy in the XY model 



Having denned the XY model and how to diagonalize it, we now proceed to find the 
entropy of the ground state as described in section 12.11 As mentioned, the ground state 
is 1 00 ■ ■ • 0) in terms of the delocalized quasi- fermions. We define the correlation matrix of 
the Majorana fermions 

T tj = (jaj -Ijli). (5.34) 

Note that the Majorana fermions have a two-to-one correspondence to the real fermions, 
so tracing out a real fermion would amount to tracing out two neighbouring Majorana 
fermions from this matrix. In the a-basis the correlator becomes 

i N 

= — (Oi t 2kOj t 2k~l — Oi t 2k~lOj t 2k) ■ 

k=l 



So far we have considered the entire system, which is in a pure state, and the entropy is 
obviously zero unless we trace out some of the system leaving us with N' < N particles in 
the subsystem. Formally, this amounts to tracing out the two columns of T corresponding 
to these particles. Now, the density matrix is assumed diagonal in some quasi-fermion basis 
denoted a. Note that this basis is different depending on the size (or shape) of the system 
traced out, but when nothing is traced out, it coincides with the basis that diagonalizes 
the Hamiltonian. Assume that in this basis the density operator can be written 

N' 

p = Y J ^n\v)(v\ Ar, = Y[[(l-n k )\ k + n k (l-\ k )} (5.35) 

r) k=l 

where rj denotes the set of fermionic occupation numbers {n\,n2, ■ ■ ■ , Wjv'} of the quasi- 
fermions, and A^ G [0, 1] are real coefficients. The sum is over all possible occupation 
numbers, = $^ m=0 " ' " Y1n n ,=o- Moreover, the ground state p = |0)(0| is identified by 
all Afc = 1. Computing the expectation value matrix f)5.34p in the a basis we find that the 
matrix can be written 

N ' ( 1 

Tij = y^ J A v {r]\[%jj]\r]) = {O i>2k O jt2k -i - O i; 2 k -iO ji2k ) f A fc - - 

r\ k=l ^ 

and in turn this means that the transformation O block- diagonalizes this matrix; 

-A fe + \ 

\ k -\ 

This shows that the eigenvalues of Y are ± (A^ — |). The idea of this exercise is that given 
the diagonalization of the Hamiltonian, the matrix V is known, and the eigenvalues can 



o T ro = i0 
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be numerically extracted efficiently even after tracing out some particles. Knowing the 
Afc's means that we know the density matrix, and can compute the entropy. The entropy 
becomes 

AT' 

S = - A, log 2 A, = H ( A *) ( 5 - 36 ) 

•q k=l 

where H(x) is the binary entropy (|2.4|) . This gives us an efficient way to compute the 
entropy of any block of real spins through the transformation into quasi-fermions. Also 
the formulas here gives an easy relation between the classical entropy of the eigenvalues of 
the density matrix in the reduced a basis and the entropy of the state itself. Figure 15.71 
shows the entropy of the XY model when half-size is traced out. The entropy clearly has 
maxima along critical lines, while the actual height of the maximum is larger along the 
7 = 0-line than the A = 1-line due to the larger central charge of this transition. 




Figure 5.7: Entropy of the XY model when half of the spins are traced out. The critical lines 
7 = and A = 1 are seen as areas of increased entanglement. Here N = 50. 
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5.4.4 Other correlators 

When we have established the matrices S and A in the Bogoliubov transformation (15.301) . 
we can compute expectation values and correlators of the model in question. Specifically, 
the second moment of annihilation and creation operators is 

N' 

(ri\a n a j n \r]) = ^ [™fc|A fcn | 2 + (1 - n k )\L kn \ 2 ] , 
k=i 

and hence follows e.g. the spin expectation value in the z-direction; 

N' 

«) = A v {2{vKal\v) ~ 1) = J2 M S fc«| 2 + (1 " ^)\A kn \ 2 } . 

rj k=l 

Some more effort gives the fourth momenta, which enables us to compute the spin-spin 
correlation function, 

This correlation will, when the distance \k — l\ is sufficiently large, fall off exponentially, 
and the correlation length £ is the cutoff of this exponential, s(k, I) ~ e~' fe ~^. Hence 
we find the correlation length in the Ising model as shown in figure 15.81 The correlation 
function follows the expected path nicely, though there obviously is no true divergence 
since the system size is finite. 
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Figure 5.8: An estimate of the correlation length in the Ising model with magnetic field. The 
phase transition at A = 1 is clearly shown as a divergence in the correlation length. The error 
bars indicate the variation of the estimate of the correlation length, thus small error bars indicate 
a good approximation to exponential decay, while the large error bars close to the critical point 
reflect the polynomial decay here. The dashed line is the thermodynamical approximation £ ~ 
|A - Ac| -1 . Data for N = 100. 
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The framework provided so far gives a good overview of the methods that can be utilized to 
find specific information about a large class of models. Nevertheless, there is an even larger 
class of models that do not conform to this framework, and where different approaches to 
finding the properties of the model has to be found. 



5.4.5 Determining criticality 

In |SO 05] we describe how to use the entanglement entropy and conformal field theory to 
determine the critical surfaces of a model. To this end, define the entropy of a block of £ 
spins in a system of size iV and inspired by the result (15. lOj) define the critical signature 
of the entropy, 

st(c) = -log 2 sin I — J . (5.37) 

This signature should fit to the actual numerical values, subtracted the value at £ = N/2 
since si/2(c) = 0. We define the error between the entropy data for a given model, §£, and 
the critical signature for a given central charge as 

t 

Thus, the natural requirement is that this error is minimal. It is obvious that the error 
would be infinite if we included the endpoints in the error, since the critical signature 
diverges when £ = or £ = N, while the actual entropy would be zero in those cases (since 
the system is empty). 

Moreover, we define an estimated central charge c es t as the central charge that minimizes 
this error, 




(5.39) 



where §£ are the measured entanglement entropies. M is the number of lattice sites that 
are within the chosen cutoff. A convenient choice is to sum over values of £ satisfying 
0.2iV <£< 0.8N, and M follows subsequently. 

It turns out that the error is very small also in near-critical systems, and a minimum 
in the error when adjusting parameters across a critical line not always corresponds very 
good to the actual phase transition. However, the if one looks at c es t when traversing a 
critical region, this almost always has a maximum at the critical point, except perhaps 
at confluence points^). Using the parameter value of that maximum and comparing to 

5 There are less well-defined critical regions, such as near the point 7 = 0, A = 1 in the XY model 
where two critical lines with different central charges intersect. Here the central charge is ill-defined, and 
the method inevitably fails. We denote these points in the parameter space confluence points. 



56 Entanglement and Its Applications in Systems with Many Degrees of Freedom 



the possible values of the central charge defined in (15.111) . one can very often determine 
the central charge and critical parameter values with great certainty. The technique is 
illustrated for a critical line in the XYZ model in Figure 15.91 
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Figure 5.9: Illustration of the technique to find a critical line as described in the text. We 
traverse the parameter line in the XYZ model with A = —1/2, 7 = 1 and A > 0. Top are si and 
the critical entropy function (|5.37|) for off-critical (A = 1, left with c es t ~ 0.368) and near-critical 
(A = 1.9, right with c es t ~ 0.512) systems for N = 10. At bottom the estimated central charge 
and the error function (|5,38p are shown when crossing the phase transition. We see the clear 
identification of the phase transition at A c = 1.9 ± 0.05. The error function's minimum is at 
a slightly higher A than the estimated central charge's maximum and the latter overshoots the 
assumed true value c = 0.5 slightly. All this is in accordance with the benchmarking done in 
[SO05], and one may conclude that the critical point is at A ~ 1.9 with the Ising central charge 
c = 1/2. 
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5.5 Criticality in the XYZ model 

We extend the XY model into what we call the XYZ model, with Hamiltonian 

H XYZ = H XY - Aj2« + v (5-40) 

n 

Despite the simple extension, this model cannot be fermionized in the way done with 
the XY model. Nevertheless, we can still find the critical regions and the corresponding 
central charges as described in Section 15.4.51 and Refs. |SO05l ISkr05b] . When scanning a 
large portion of the parameter space, using N = 10 (with symmetries) is useful, since the 
computation is rather reliable and the time consumption is reasonable. Near confluence 
points one can use larger systems for small areas of the space. All data in the next section 
are taken with N = 10 to demonstrate the power of the method using such a small system. 



5.5.1 Critical surfaces 

We investigate the XYZ model in the range — 2 < A < 1/2, A > 0, and 7 > 0. First, 
consider the case 7 = 0, which is critical in the region shown in Figure 15.101 That is, in the 
region where A < 1 — 2A and A larger than a lower limiting line as indicated in the figure. 




Figure 5.10: The critical region of the XYZ model with 7 = 0. The critical surface has central 
charge c = 1, and is faded where the critical nature of the model becomes indeterminate for this 
determination which is done with N = 10. The black line indicates the intersection with the XY 
model. The model is symmetric under A — ► —A. 
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In the region A « -1 and |A| < 0.9 the critical region ends in the sense that for smaller 
A the model is massive, but the end of the criticality is indeterminate. This means that 
the estimated central charge falls below one, and the minimizing error gradually becomes 
larger, though there is no sharp transition to distinguish the critical from the non-critical 
region. The two lines mentioned above, however, have very sharp transitions between 
critical and non-critical regions and the central charge is readily identified at c = 1. The 
c = 1 critical surface exclusively belongs to the regime 7 = 0, except perhaps for A < — 1, 
where this surface seems to spread out into regions of non-zero 7. The details here are not 
investigated in detail, though. 

In addition to the c = 1 surface, there is a c = 1/2 surface in the parameter space, which 
in the intersection with the XY model at A = is at |A| = 1 for all 7s. When A is 
larger than 1/2, the c = 1 critical surface vanish for 7 = 0, while the c = 1/2-line becomes 
indeterminate. However for A < 1/2, the c = 1/2 transition is illustrated in Figure 13. Ill 
It has a very weak dependency on 7, and only for positive A. 
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Figure 5.11: Intersections of the critical surface with c = 1/2 in the XY Z model for three values 
of 7 seen in the A- A plane. There is an uncertainty in the A value at each point of about ±10 -2 
due to the inherent uncertainty in the method. The lines are computed based on the maximal 
estimated central charge for iV = 10. The horizontal line is the intersection with the XY model. 
Note that it is possible to draw the line for 7 = 0.5 only for A < 1, since it becomes indeterminate 
for this system size for larger A. This apparent indeterminacy is due to the proximity to the c = 1 
critical line, and larger system sizes will make it possible to draw this line further. The thick line 
is A = 1 — 2A, the end point of the c = 1 surface, and presumed starting point for the c = 1/2 
surface for 7 — ► 0. 



Quantum Critical Systems 59 



5.6 Time evolution 



Given the formalism of the fermionization developed in section I5.4[ it is easy to find the 
time evolution of the state [Skr06b]. In the Heisenberg picture, the Majorana fermions 
defined in (I5.26p has a time evolution according to 

-^7 fc = i[H, %} = -2i (5-41) 

i 

This is a linear first order system of differential equations, whose solution is 

lk{t) = T kl % (5.42) 
where jk — %(t — 0). The time evolution matrix T is given by 

T H (t) = ^S fc ,S*e- 2i ^ (5.43) 

i 

when S is the eigenvector matrix of C and £ are the eigenvalues. 

A state that is in a thermal equilibrium with its environment at an inverse temperature (3 
is a mixed state described by its density matrix, 

a = l e ~^ (5.44) 

where H is the system's Hamiltonian and Z = Tre~@ H is the partition function. Given a 
basis in which the Hamiltonian is diagonal H = Ylk -^kl^k) {ipk\, this becomes 

<x=|^e-^|V>*>(V>*i. 
k 



In |Skr06bj we compute the time evolution of an excited pure state in an Ising chain with 
a magnetic impurity, and show that the eigenvalues of the reduced density matrix for a 
part of the chain resembles those for a state in thermal equilibrium at some temperature 
To this end, we compute the classical fidelity between the eigenvalues. If we denote 
the ordered eigenvalues of the reduced density matrix and the thermal density matrix 
ip k = e~@ Ek /Z, the classical fidelity becomes 

f c (x, ip) = 22 V^kPk- 

n 

The classical fidelity is an adequate measure of the distance between two probability dis- 
tributions, being 1 if the distributions are equal, and less than one else. Hence the fidelity 
is not a metric, but is intimately related to the trace distance, which is, and thus can be 
considered equivalent to a metric |NC00] . 
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For two arbitrary quantum states described by density matrices p and a, the fidelity 
between them is defined as 

F(p, a) = Tr y/p^ap 1 / 2 , (5.45) 

which reduces to the classical fidelity of their eigenvalues in the case where p and a com- 
mute. However, for the non-commuting case, it gets more complicated. Given that the 
two states can be fermionized, and thus written in terms of the Majorana fermions 

/,-/"/,-"/ ! N p - : 




n*(i+e- n *o 

1 



we need to compute the matrix product The Campbell-Baker-Hausdorff formula 

for exponentials of non-commuting variables A and B reads 

e^e 8 = e v where 

V = A + B + ±[A,B\ + ±[A-B,[AM + ---, (5-46) 

given that the series converges (which is not obvious). If we have operators = 2%7iTj 
and 5 = |^7i7j f° r some anti-symmetric matrices 9 and £, we can use the commutator- 
anti-commutator relation 

to find that 

[6,2] =-<p ij n /i y j =$ (5.47) 

where <p = [0,£\. This relation holds for all iterated commutators involved in f )5.46p . and 
we may therefore conclude that 

e A e B = e v , 

where A = ^Aijjfij, B = ^Bijjijj, and V = ^ijlilj with D given by the expansion 
( 15.46!) . The generalization to more than two factors is simple, and we can compute the 
matrix product a l / 2 pa 1 ^ 2 and (block) diagonalize the result into the matrix exp (— i^i/TiT?) 
with new Majorana fermions 7. Thus the fidelity becomes 

F{a,p)=Tcy/Njr p e-^ Ki ^, 

which is straightforward to compute. 

Hence it is possible to find the full fidelity rather than just the classical fidelity of the 
eigenvalues, with the possible catch that the Campbell-Baker-Haudorff series (I5.46P may 
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not converge. However, it is not expected that the difference between the classical and 
quantum versions will be large, since the diagonalization of the Hamiltonian and the den- 
sity matrix will be "almost" the same. That is, the diagonalization matrices O a and O p 
respectively fulfill 0~ 1 O p ~ 1. However, discrepancies in this equality may map back to 
small eigenvalues that are of little importance to the fidelity. 



5.7 The quantum Hall effects 

The exposure of an essentially two dimensional conductor to an magnetic field gives rise 
to a voltage perpendicular to the current flowing through the slab, which has long been 
recognized as the Hall effect. A schematic experimental setup is shown in Figure 15.121 The 
Hall resistance R xy is the Hall voltage across the slab V y by the exposed current I x , and in 
the classical case the Hall voltage is proportional to the imposed magnetic field. Indeed, 
straightforward calculation gives the law R xy = B/ne, B being the perpendicular magnetic 
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B [T] 

Figure 5.12: The resistivity of a quantum Hall system, sketched in the upper left. Both resistiv- 
ities in orthogonal R xx = V x /I x and transverse R xy = V y /I x are shown and the most important 
fractional and integer quantum Hall states are indicated. The experiment uses a 2D system with 
electron density n = 2.33 x 10 15 m -2 cooled to 85mK. Image from |STG99j . 
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field, n the density of conduction electrons^ and e the elementary charge. When the system 
is cooled down, and the magnetic field becomes large, the Hall resistance becomes quantized 
at exact levels of 

_NJi_N 1 _ eB 

V eB N s ' s h ' 

with v an integer called the filling factor, as first reported in 1980 by von Klitzing et al. 
[KDP80]. h is Planck's constant, and the phrase filling factor refers to the fact that N s 
is the number of available states in each Landau level and N e the number of available 
electrons. As the magnetic field is increased, N s increases proportionally, and when v 
becomes integer, the Fermi level resides between Landau levels, thus all available states 
are occupied and the electron fluid is essentially incompressible until the magnetic field 
increases to another integer. This integer quantum Hall effect (IQHE) is explained as a 
problem of non-interacting (apart from the Fermi statistics) electron problem first explained 
by Laughlin in 1981 |Lau81j . 

In 1982 Tsui et al. |TSG82] identified levels of constant resistance also at levels of fractional 
z/, a result completely unexpected at the time, and a matter much harder to explain than 
IQHE. The fractional quantum Hall effect (FQHE) occurs at filling factors v = p/q for 
certain fractions, where the most prominent levels are those with v = 1/q and v = 1 — 1/q 
with q odd (called primary states), and v = p/(2p ± 1) (higher order states), cf. Figure 
15.121 However, many more fractions appear, even fractions with even denominator. FQHE 
is, in contrast to the integer variant a complex many-particle problem, and a number of 
explanations exist. Perhaps the most prominent of these are Laughlin's trial wave function 
[Lau83j . which for the v = 1/q state of n particles reads 



n / 1 n 

Zn) = II (** ~ Z ^ 6XP ~i/2 Yl 

i<j \ k=l 



\ z k\ 



where Zk is the 2D complex position of particle k, and normalization is omitted. Iq = 
a/ h/eB is the magnetic length. This reflects attaching q flux quanta to each electron, 
since the charge concentration of electrons and charge deficit of a vortex associated with a 
flux quantum attract each other, thus creating fractionally charged excitations. 

More recently, Jain's composite fermion picture [Jai89j has become a more fashionable 
explanation, particularly since it also explains many higher-order states. Essentially, the 
system of strongly interacting electrons in a strong magnetic field is mapped onto a system 
with a reduced magnetic field and a series of non-interacting composite fermions consisting 
of one electron and two flux quanta |Jai00j . Thus, the fractional quantum Hall effect 
becomes a integer quantum Hall problem of fractionally charged fermions. 



6 In a typical setup on uses GaAs/AlGaAs heterostructures with n ~ 10 15 m 2 . [STG99] 
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5.7.1 FQHE on a torus 



Most of the literature has so far focused on explaining the inner workings of the FQHE, 
and little is known about the entanglement properties of the quantum Hall states. Nev- 
ertheless, it has been argued that on the IQHE plateaus the entanglement is zero, since 
this is essentially a non-interacting system [Shi04j. As the magnetic field is increased (or 
decreased) and the plateau changes from one FQHE state to another, we have a quantum 
phase transition, and possibly the criticality technique described in |SO0 5j and Section 
15.4.51 could be applied here. 

Even though the QHE system is a generic two dimensional problem, it can be mapped 
onto a ID problem [CP95, YHL83J. That is, with periodic boundary conditions, one have 
N different single-particle wave functions in each Landau level, which in a rectangular 
geometry with sides L x and L y are 



1/2 00 r. 1 

exp -2 (Xj + kL x )y - -^(Xj + kL x - x) 2 



(5.48) 



Here, 1 < j < N s and Xj = jL x /N s are localizations along the a;-directioE0. This amounts 
to N s localizations along the x direction with Gaussian with of the order of the magnetic 
length. The probability distribution is independent of y. 

In a second quantized version, the Hamiltonian is 

H = J2 Wa)aj + Y Y Y A hj2j3j4<4i a h a h a j4 (5-49) 
j ji h is 34 

where dj is the destruction operator at site j. The single electron energy W is a known 
constant |YHL83] and the coupling term is 



ne 2 



e L x L y ^ q q 



exp 



-~llq 2 - \q x L x (jx -j 3 )/N s 



r/ r/ 

ji ~34 ,q v L y /27T ji + j 2 ,jz +j4 



where 5' is a Kronecker delta modulus N s and V(r) is the real-space Coulomb interaction. 
The sum over q is a sum over all allowed q- vectors (except the zero vector), that is with 

integers s and t such that q = ^^p, x^J- Thus a basis of the Hilbert space is defined by 

the A^ occupation numbers, \ji, - ■ ■ , Jjv )- For a given filling fraction v this means that the 
Hilbert space has dimension (^ s ) , which obviously quickly becomes far too large to handle. 
Indeed, the Hilbert space's dimension grows exponentially with the number of electrons for 
a fixed is, and faster than 2 e . Again, utilizing symmetries makes it possible to go slightly 



7 This example is in the Landau gauge A — Bxy. 
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further. With the most naive approach however, investigating six particles in a v = 1/3 
system involves matrices of size (g 8 ) = 18564, which is computationally time consuming. 
Larger systems are essentially inaccessible. 



5.7.2 Entanglement in the FQHE 



The entanglement in the FQHE case is conceptually different from the spin chains setting. 
This is mainly due to the complications involved in dealing with identical particles. In 
spin chains we consider each particle fixed with a non-overlapping wave function, and the 
particles can hence be thought of as separate, or non-identical. This is not the case in the 
FQHE where the overlap of the wave function may be considerable. Entanglement is ill- 
defined in this case, mainly since the anti-symmetrized wave function defined by the Slater 
determinant already carries quantum correlations while they particles involved cannot be 
said to be entangled, in the sense that there exists no less correlated state. It has been 
suggested that entropy of reduced density matrices larger than that of a Slater determinant 
state is an entangled fermion state [GM04| . For bosons the argument is likewise, but 
since the least correlated bosonic state is a product state whose entropy is zero, positive 
entropy indicates an entangled state here. However, it is not obvious that the entropy 
(possibly subtracted the Slater determinant entropy) is a genuine entanglement measure. 
The entropy of a single particle in a Slater determinant with N particles is log 2 N. 

For the FQHE case the reduced density matrix of N' e < N e electrons is of size (jj) and thus 
the reduced density matrix may be larger than the original matrix, if N' e > N e /2. Despite 
this apparent paradox, it is possible to compute the entropy of the reduced density matrix. 
To this end, consider the density matrix of the system in the basis {|110), 1 101) , |011)}, 
P = iPij] m this basis. Tracing out one particle to remain with one amounts to tracing over 
all wave functions with one particle, and in the new basis {|100), 1 010) , |001)} this becomes 



(l..|p|l..) + (.l.|p|.l-) + (--l|p| 



Pll + P22 
P32 
P31 



P23 
Pll + P33 
P21 



Pl3 
Pl2 
p22 + P33 



Numerical calculation of the entropy for the ground state of (15.491) by the prescription 
above, shows that the integer quantum Hall effect has entropy S n = log 2 (^ s ) when the 
system is traced out to leave n particles, as predicted in [ILO06] . Tracing out to a single 
particle, the entropy is simply S\ = \og 2 N e . The values for the entropy in the IQHE 
states are independent of the nature of the interaction, since IQHE is an essentially non- 
interacting phenomenon. 

For the FQHE the case is more complicated, and exact analytical results are rare. Using 
the prescription above, some basic results for very few particles can be found. In particular, 
we find that in the limit where L y <C L x , the entropy is exactly that of the IQHE, log 2 (^ s ) 
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when n particles are traced out. In this limit, the localizations are so far apart that the 
particles are essentially non-interacting, and the IQHE regime is recovered. However, as 
the aspect ratio L y /L x is increased, at some value the entropy for a single particle jumps 
discontinuously to a higher value, and after this increases slowly with the aspect ratio until 
saturating. 
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